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ELLIPTIC FUNCTIONS AND INTEOAIS WISH REAL 
MODULUS IN FLUID MECHANICS* 

By Robert Legendre 

SUMMARY 


Advantage of the elliptic functions and of the more general functions 
of Schwarz for fluid mechanics. Flows outside and inside polygons. 
Application to the calculation of an elbow diffuser for a wind tunnel. 
Properties of the elliptic integrals of the first kind and of the elliptic 
functions. Properties of the theta functions and decomposition of the 
elliptic functions into products of theta functions. Properties of the 
zeta functions. Decomposition of the elliptic functions into sums of 
zeta functions and calculation of the elliptic integrals. Applications 
to the calculation of wing profiles, of compressor profiles, and to the 
study of the vibrations of airplane wings and of compressor vanes. 

The manuscript of the present paper was checked by Mr. Eichelbrenner 
who corrected several imperfections and suggested numerous improvements 
to make reading of the paper easier. However, 'the limited subject does 
not permit filling In more than an incomplete knowledge of the properties 
of analytic functions. 


INTRODUCTION 


The solutions of a very large -number of problems in fluid mechanics 
are expressed with the aid of elliptic functions. The mechanism of the 
role of these functions is not difficult to analyze. 

From one point of view, the elliptic functions can be considered as 
the simplest ones (after the exponential and circular ones which they 
generalize) among the solutions of differential equations whose coeffi- 
cients are polynomials. Thus, it is natural that one must resort to 
elliptic functions when a somewhat close approximation is desired. 

From another more geometrical view point, the majority of the prob- 
lems which can be solved with the aid of exponential and circular func- 
tions are related to schemes which assemble the given quantities on a 

*"Les Fonctions et Integrales Elliptiques a Module Reel en Mecanique 
des Fluides," ONERA, Publication No. 71, 195*K 

NACA Reviewer’s note: The original French publication contains cer- 

tain typographical errors and obvious omissions in equations that have 
been corrected without comment. 
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unique segment or else on a unique curve or on a family of curves, deduced 
one from another "by a simple periodicity. As soon as the given parameters 
are on two separate curves or on a family of curves deduced one from 
another by a double periodicity, the elliptic functions are introduced. 

The elliptic functions are analytic functions, the field of which 
represents directly the plane flow of an ideal incompressible fluid. 

The conformal transformations allow the association of the fields of flow 
around very different obstacles . On the other hand, the study of the 
flow of compressible fluids is largely based upon the acquired knowledge 
of analytic functions, at least as far as search for approximations by 
various artifices is concerned. 

Among these artifices one must mention the study of the hodograph, 
that is, of the potential of the plane flow of a compressible fluid repre- 
sented with the aid of the velocity. The relation between the potential 
and the velocity is in fact much closer to the one expressed by an ana- 
lytic function than the relation between the potential and the coordinates 
in the physical plane. The theory of conical flows satisfying an equation 
which is linearized by approximation, on the other hand, leads to the 
study of the three-dimensional flow of a compressible fluid in terms of 
the plane flow of an incompressible fluid. Finally, the projection of 
the velocity of the almost uniform flow of a compressible fluid onto a 
plane perpendicular to the mean direction of this velocity is approxi- 
mately the velocity of the flow of an incompressible fluid. The nearer 
the Mach number of the mean flow is to unity, the closer is the 
approximation. 

The theory of the elliptic functions is generally very little known 
among engineers, not so much for want of mathematical treatises as because 
of the lack of work establishing the connection between these functions 
and the flows which they can represent. One finds quite often a thoroughly 
documented treatise which does not contain a single figure. 

The present paper uses largely the geometrical methods taught in 
1930 at the Ecole d’Application du Genie Maritime by the Institute mem- 
ber, Mr. Emile Barrillon. It does not attempt the rigorousness to be found 
in the mathematical treatments but the intuitive and quick solution of 
the problems . The general theorems relative to the analytic functions 
permit a justification of the exactness and uniqueness of the solutions 
found. 

Although the paper belongs in the domain of applied mathematics, it 
establishes some properties of the elliptic- functions which are, to our 
knowledge, original. 
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1. SCHWARZ* TRANSFORMATION 


1.0 Flow Outside or Inside a Polygon 


The later developments will show that the elliptic functions with 
real modiolus and the integrals -which can "be attached to them represent 
flows outside and inside quadrirectangular polygons. They generalize 
the circular functions which are associated with flows limited by 
birectangular polygons and are themselves particular cases of the Schwarz 
functions defining the flows outside or inside arbitrary polygons. 

Finally, Schwarz* functions are the simplest of the automorphic (or 
Fuchsian) functions of Poincar£. 

Without aspiring to such large generalizations for the limited sub- 
ject of the present paper, although the geometrical methods are of a 
character to facilitate the understanding and to generalize the applica- 
tion of Poincare's functions, it is useful to associate the elliptic 
functions with real modulus to the functions of Schwarz. 

Let us recall the principal properties of the analytic functions 
which will be useful later on. 

An analytic function is associated with a representative field 
constituted by a net of two families of orthogonal curves forming a 
grid. Each curve corresponds to a constant value of the real part or 
of the imaginary part of the function. For the interpretation in terms 
of fluid mechanics, one of the families corresponds to the equipotentials, 
and the other to the stream lines. 

An analytic function of a variable, which Itself is an analytic 
function of a second variable, is an analytic function of the second 
variable. Its representative field in the plane of the second variable 
may be constructed by a point-by-point conformal correspondence, that 
is, locally conserving the angles and the ratios of the lengths except 
of the second order. An analytic function such as the one above which 
defines the correspondence between the two variables is thus associated 
with a conformal geometrical transformation. 

An analytic function defined in the entire plane of a variable dif- 
fers only by a constant from the sum of its principal parts in the 
neighborhood of all its singular points. If an analytic function is 
defined in a limited region of the plane, but if It is real or has a con- 
stant imaginary part on the contour of this region, and if it is possible 
to make this region correspond, by conformal transformation, to a half 
plane or to the interior of a circle in such a manner that the function 
can be defined by analytic continuation in the entire plane of the new 
variable, it differs only by a constant from the sum of its principal 



k 


NACA TM 


parts in the neighborhood of all its singularities and from the singu- 
larities of the analytic continuation in the entire plane of the new 
variable. . 


1.1 Hodograph of the Flow Outside or Inside a Polygon 

In accordance with the foregoing review, it suffices to characterize 
one flow defining a conformal transformation of the polygon into a circle 
or into a straight line in order to determine all flows containing doublets, 
sources, sinks, or vortices In limited number or distributed on the curves 
In Infinite number. These flows will, in fact, be associated with an 
analytic function, determined except for one constant for these enumerated 
singularities and their images, by Inversion with respect to the circle 
or by symmetry with respect to the straight line. 

For reasons of symmetry, the retained fundamental flow will be that 
of a circulation around the polygon or that of an isolated vortex Inside 
the polygon. _ 

Let u be the complex variable in the plane of the polygon and 
x(u) the complex potential of the flow whose real part Is the potential 
of the velocities, while a streamline corresponds to & constant value of 
the imaginary part. We can select X real on the contour of the polygon 
(figs. 1 and 2). 

The velocity has a constant direction on every side of the polygon. 

This property may be easily characterized by use of the hodograph, that 
is to say, of the velocity or also of Its Inverse, studied as a function 
of X - 


I - pe ie ■ too 

where p is the Inverse of the intensity of the velocity, and 0 the 
angle of Inclination of the velocity, constant on one side of the polygon. 

Even more convenient will be the study of the logarithmic hodograph 


ln£ = lnp + 10 

the imaginary part of which is a constant on one side of the polygon cor- 
responding, in addition, to the real axis of the X- plane or to the circle 

iX 

of radius 1 in the plane of the variable z = e . The-variable z will 
always be associated with the function X In the following calculations . 
Comparison of the figures 1 or 2 with figure 1(a) shows, in fact, that the 
con formal transformation u(z) defined implicitly for X(u) and X(z) 
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makes the exterior or interior of a polygon correspond to the interior 
of a circle. Another transformation u(t) , defined further on, vill make 
the exterior or interior of a polygon correspond to a half plane . 


1.2 General Solution 

Let us choose to study Inf as a function of z = e iX . The vertices 
of the polygon correspond to the unknown points of the unit circle 
zi,z 2 . . ., 2 ^. • which are singular points of the function u(z). 

No other singular point exists, except perhaps at the origin, center of 
the circle, which we agree to make correspond to the point at infinity 
of the field outside the polygon or to the point of reference of the 
field inside the polygon since the desired circulation has no other 
singular point outside the polygon but the point at infinity in the case 
of external flow, and the center of the chosen vortex in the case of 
internal flow. 


The constancy of the imaginary part of Inf on the unit circle in 
the plane of the variable z permits defining the analytic continuation 
of the function toward the outside of the circle. It is sufficient to 
make correspond to z and Inf, for an internal point, the relation 
between l/z and Inf + C"k e defining an analytic function at the exter- 
nal point which is the inverse of the first with respect to the circle. 
The two f inactions join on the circle where z = l/z and 

Inf = lnp - 10 

It suffices to adopt as constant 210, that is to say, the angle of 
inclination of the side of the polygon multiplied by 2i, arranging cuts 
between the singular points and, for instance, the small-straight lines 
extending the radii (fig. 3). 


The mode of analytic continuation shows that there exists no singu- 
lar point outside the circle, except perhaps at infinity. All singular 
points are thus known and the function may be defined by the sum of its 
singularities . 


In the proximity of a vertex u^. of the polygon, the function X(u) 
behaves like 


X = + C-^u - u^j* + . . . 


where 0^ is the angle at the vertex measured toward the fluid. 
Consequently, 
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u - u k 


; (x - \) 


s - g-°3(x-x k )” 


0 , 

Jt- 1 


Since z(X) = e iX is regular at a vertex z k = e iXk , 

i-i . 

S = Ci f (z - z k ) n + . . . 


lag 


= - lj log (z - ZJt) + . . . 


At infinity, for the flow outside a polygon, and around the selected 
vortex for the flow inside a polygon, the function X(u) behaves like a 
logarithm 


X = ±1 in u + . . . 

The plus sign corresponds to a flow outside a polygon and the minus 
sign to a flow inside one. 


1 _ dX 
g du 



tie 


TiX 


In g = tiX + . . . = +ln z 


It is not necessary to study the infinity of the z -plane where the 
singularity is determined by the analytic continuation 

Finally, lng is defined by the sum of its singularities 

In g = tin z + Z ^ - lj In (z - z k ) + C te 

If one notes that the sum of the deviations of the- stream on the 
contour of the polygon is four right angles, then 


- 1 ) - * 2 
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where the plus sign corresponds to the flow outside a polygon. Con- 
sequently, the function ln£ behaves, for an infinite z, like ±lnz. 
On the other hand 


z 





and 



The flow outside or inside a regular polygon with n sides corre- 
sponds to real values of X stepped at 2k / n and to values of 

exactly equal to rf ± 2it/n. The plus sign corresponds to the external 
flow. 

1.3 Extension to More Complex Polygons 

The formula established above for values of 0^ between 0 and 2k 
can be generalized. 

A zero value of corresponds to a polygon vertex at infinity 

between two parallel sides (fig. 5)» 

A value of 0^ equal to 2k corresponds to a point of return 
(fig. 6). 

A negative value of corresponds to two infinite branches 

inclined by -0^. (fig. 7)* 

A value of 0^ higher than 2k corresponds to an overlapping 
between the two adjoining sides (fig. 8). 
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Finally, the formula defines a polygon, for a real X , if 
is a positive or negative integer different from *1. 


- 1 ) 


Even when all values of % are contained between 0 and 2 rt, the 
choice of the values of for the representation of an a priori given 

polygon is difficult. It is generally expedient to resort to the elec- 
trical analogy, and. the analytic expression is no longer applicable except 
to the exact numerical study in the neighborhood of the singularities . 


When X varies from 0 to 2rmt where m is the entire value of 

i z(— - l) , the variable u resumes its initial value and the polygon 

closes. In fact, all cuts in the plane of the variable z may be chosen 
outside the unit circle, and the integral of £ dX on that circle is equal 

to the integral around an infinitely small circle. Since £ is, except 
for one factor, equivalent to z, the integral 


u = 



is itself equivalent to z, except for one factor, and, since it is inde- 
pendent of the integration contour, it is necessarily zero. 


1.4 Change of Reference 

For the flow inside a polygon and for the flow limited by a polygon 
having infinite branches, the vortex of reference was chosen arbitrarily. 

If another field or reference X’ is selected, it may be useful to 
determine its relation with the initial reference X. 

In the plane of the variable z = e* - *, the field X* is that of a 
vortex inside the unit circle. It is, except for one constant, defined 
by this singularity situated at a point the complex variable of which 

will be denoted by and by the image of the latter with respect 

to the circle. 

It is convenient to make a change of axes which makes the symmetry 
evident, taking ze" 1 ^ as the variable. The new function X* then is 

X' = -i In (ze - ^ - e” a ) + i ln^ze ^ + 0' - ia 
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The imaginary part of the constant -ia was chosen so that X* 
should be real on the unit circle in the plane of z, image of the con- 
tour of the polygon. 

i y * 

If one sets z' = e , the relation "between z' and z is 
tomographic 

(z 1 e -1 ^') (ze“ iP ) - e a ^z , e -i ^' + ze” ip ) + 1 = 0 

The relation between X 1 and X may be put in another form, con- 
venient for the real values 

, X T - p f , X - p , , , a ^ 

tan — tan — + tanh — = 0 

2 2 2 

This expression shows the advantage of the function 

t(u,p) = tan - - = -i z - ~ — 

2 z + 

which is real on the contour of the polygon and depends on a real param- 
eter p. The above formula shows that t(z,p) is represented in the 
z-plane by a doublet on the circumference at the point z = -e^-0, assuming 
the unit circle as a streamline . The conformal transformation z(u) thus 
makes the function t(u) correspond to t(z)j this function t(u) is 
represented in the field of the variable u by a doublet on a side or at 
a vertex of the polygon (fig* 10), according to whether the value chosen 
for p does not or does correspond to the vertex of the polygon in the 
transformation z(u) . 

When no reason of symmetry makes it advisable to prefer the func- 
i X 

tions X or z = e for defining the conformal transformation of a 
polygon into a circle, it is often convenient to utilize a transforma- 
tion t(u,p) for a judicious value of p which makes the real axis of 
t correspond to the polygon, and one of the half -planes limited by the 
real axis to the interior or to the exterior of the polygon. The trans- 
formation t(u,p) may be defined directly by 



where 
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This expression Is suitable only when the doublet of t(z) has not 
been chosen at a vertex of the polygon. If, on the contrary, 0 *= Xj + it, 

the expression maintains the same form, but the term corresponding to ' 
must be omitted in each of the products. 

a/2 

The change of reference t = tanh ^ — is .particularly convenient 

with the latter expressions. It shows that has the same form as 

d dt ' _ . 

— : which was obviously necessary, 
dt 

For the flow inside a polygon, £ - lj = -2 and the term in 

p 

1 + t disappears. In the other cases, the presence of this term indi- 
cates the existence of a point at Infinity or of a critical point at the 
center of the vortex X(u). 


If all values of 0j£ are Integer multiples of TI, the expression 

is a rational fraction of t and the integral Is expressed with the 

aid of rational fractions and of logarithms, that is, of functions which 
are the inverse of exponential or circular functions. 

If one of the values of 0^ is an odd multiple of rc/2, there exists 
at least one other value of of this type If all others are multiples 

of H, since the sum of the values- of 0^ Is a multiple of 2rt. The 

polygon has then two right angles and du/dt which is a square root of 
a polynomial of the second degree, multiplied by a rational fraction in 
t, Is again integrated with the aid of rational functions and of exponen- 
tial or circular functions. 

If three values of are odd multiples of rt/2, there exists a 

fourth value of 0^ of this type for all others to be jnultiples of jr. 

The polygon has then four right angles, and 'du/dt Is the square root of 
a polynomial of the fourth degree In t, multiplied by a rational function. 
The function u(t) Is called an elliptic integral and its properties will 
be studied in the following chapters. ~ 

1.5 Field of Doublets or of Sources and Vortices in a Polygon 

With the use of a. conformal transformation which makes a fraction of 
the plane limited by a polygon correspond to the interior of a circle or 
to a half plane, it is easy to study the field F(u) of doublets or of 
sources and sinks limited by a polygon which is transformed into a field 
defined in the entire plane by analytic continuation with the aid of images 
with respect to the circle or with respect to the real axis. 
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We shall choose to utilize a transformat Ion t(u) for studying the 
field F(u) -which is transformed into F(t) defined In the entire plane 
by symmetry around the real axis of t. 

The Slim of the singularities of doublets, or of poles of a higher 
order forms a rational fraction F^(t ) the analytic continuation of 

which is Fi(t), obtained by replacing the coefficients by their con- 
jugates. The field F(t) = Fj_(t) + Fj_(t) is a rational fraction with 

real coefficients. Inversely, such a rational fraction defines by a 
transformation t(u) a field of poles F(u) admitting a polygonal 
streamline . 

The singularities corresponding to sources, sinks, and vortices are 
logarithms. The function F(u) becomes, after the transformation t(u) 

F(t) = z[A ln(t - tj) + A ln(t - tj)J + C te 

Inversely, an expression of this form defines a field of vortex 
sources which contains, except when Z(A + A) = 0, a source on the con- 
tour of the polygon, for infinite t. 

A frequently encountered particular case corresponds to the sources 
and sinks, the strengths of which are equal or are in a simple fractional 
relationship. The function F(t) is then the logarithm of a rational 
fraction. 

The field of a line of sources or of vortices is defined by 


F(t) = J' A(s)ln(t - t s )ds + J' A(s)ln(t - i 8 )di 


where A(s) is an arbitrary intensity distributed on the line C as a 
function of a parameter s which can be the curvilinear abscissa. The 
intensity A(s)ds is real for sources and purely imaginary for vortices. 


1.6 Analytic Continuation, Periods, Case of Reduction 

The analytic continuation of a real function (or one with a constant 
imaginary part) on each of the sides of a polygon is obtained by successive 
symmetries with respect to the sides. The continuation thus defined is 
generally multiform and the procedure followed depends on the order of the 
symmetries. 

The continuation is periodical. Two successive symmetries, one with 
respect to a side of the polygon, the other with respect to the homologue 



12 


MCA TM 1^35 


of another side in the first symmetry, are equivalent to the combination 
of a translation and of a rotation by an angle double that of the two 
sides . 

If the symmetries are combined in such a manner that the rotations 
have as their sum a multiple of 2rt without the translation being zero, 
the function F(u) admits this translation as the period. 

Figure 11 illustrates the construction of a period in the case of a 
triangle. 


The rotation after the sixth symmetry is 


-2a - 2p - 27 = -2n: 

For the rectangles and the equilateral, triangles, the symmetries 
show the way, and the sense of description which is inverted after each 
symmetry, is reproduced after rotation around a vertex. A function F(u), 
which is real on the contour of the rectangle or of the triangle, is then 
uniformly continued. For regular hexagons, -the symmetries again show the 
way but two turns must be described around a vertex in order to reproduce 
the sense of description. We shall find finally that the uniformity of 
the function defined in the entire plane, combined with the double perio- 
dicity, is characteristic of elliptic functions. 

The analytic continuation by symmetries permits utilizing, for the 
study of a symmetrical flow in a polygon, the transformation relative to 
the half polygon limited by the axis of symmetry. In order to study, for 
instance, the field F(u) of a doublet on the axis of symmetry of a 
quadrilateral with a circulation such that the stagnation points are at 
the vertices (fig. 12(a)), it will be convenient to utilize the trans- 
formation t(u) relative to one of the triangles (fig. 12(b)) and con- 
tinued analytically in the other 

£_1 Ss - 1 

jjf = (t - if (t + l) rt 

In the plane of the variable t, the function F(t) has the appear- 
ance represented by figure 12(c). The cuts will be eliminated by the 
classical transformation 2t = s + l/s and .the function F(s), uniform 
in the entire plane, will be defined by its singularities 


(s) = ■ . + is infs - e i<0 ) + - A g . " 1 ^ . 13 ln(s - e” ia) ) + C 

s - e 1 " V ' e _ V 


s - e 


The condition of separation of the stream at C in the F(u)-plane 
imposes an infinite branch in the F(t)-plahe "and a stagnation point at 
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the origin in the P(s) -plane, either with dF/ds being zero for s = 0 
or else B = -A cot a>. The constants oi, A, C, remain arbitrary. The 
first defines the position of the doublet in the F(u)-plane, the others 
fix the scale and the origin of F. Finally, the field F(u) to be 
studied is defined parametrically by F(s) calculated above and u(t) 
defined by integration of du/dt whose variable t is itself a function 
of s. 


Another particular case of reduction pertains to the flows limited 
by an equilateral triangle which can be related to flows inside a particu- 
lar rectangle, that is, as we shall see, to the elliptic functions and 
integrals . 


In fact, we shall study the field of a doublet at one of the vertices 
of the triangle, corresponding to a function t(u) defined by 



If we set 


then 


1 


- t2 = 



and y(u) defines , as the later chapters will show, an elliptic integral 
with a real modulus. 


In general, numerous cases of reduction of the functions of Schwarz 
exist which may be found by analytical or geometrical methods. The sub- 
ject of the present report is limited to the study of the reductions of 
the elliptic functions and integrals with real modulus which will form 
the subject of the following chapters. 


1.7 Examples of Application 

Assume that turning vanes for a right-angle elbow diffuser intended 
for a wind tunnel have to be defined. 

One intends to construct the vanes very simply by means of one metal 
sheet reinforced with the aid of a second, soldered to the first (fig. 13 )* 
For the setup of the calculation, the thickness of the sheets is neglected 
and the regions of the trailing edge and the leading edge are assumed to be 
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straight. It is desired, moreover, to orient the leading edge so as to 
avoid its becoming bent. 

Let F(z) be the complex potential of the flows (fig. 14) . To 
choose the distribution of velocities, it is convenient to construct the 
inverse hodograph 

t (f ) = — 

* dF 

Let us follow the streamline -which divides at the leading edge 
(figs. 14 and 15). 

From B to C, the direction is constant and £ varies on a radius GBC 
inclined by e. 

From C to D, the direction of (; changes and its intensity increases. 
Its extremity describes the curve CD. 

From D to E, the direction is again constant and £ describes the 
radius DE, inclined by + e 1 . 

The same reasoning applied to the upper surface of the vane fixes 
as limits of the field F(£) a curvilinear quadrilateral CDD'C'. The 
variation of the velocity will be quite continuous if logarithmic spirals 
are chosen for the curves CD and C'D*. 

Let us now carry out the transformation ln£. The curvilinear 
quadrilateral is transformed into a parallelogram (fig. 16 (a)), which 
itself corresponds to a circle or to a straight line by means of a Schwarz 
transformation (fig. 16 (b)) 

= C 0 (sin X) 2 + 5 sin(Xj_ - X) 2 

where the real part is chosen for 0 < X < X^ and where C Q is a real 

ponstant determined by the condition that the Imaginary variation of ln£ 
between C and D is 

i(| + e + e') 

pjt -i+5^. _ « i “ 8 

^ + e + e' = C Q cos 5 J [sin xj |sin(x-j_ - XjJ &X 
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In order to choose the value of ^ corresponding to the image of 

the point at infinity upstream., one will make X vary hy imaginary values, 
starting from an arbitrary real value Xg such that 0 < Xg < X^ 


ie- = 




stopping the integration at the value X = Xg + irjg such that this equa- 
tion is satisfied. One will operate in the same manner for defining Xq., 

corresponding to the image at infinity downstream, starting from a 
value rt + Xj such that 0 < X^ < X^ 


ie’ 





dX 


Practically, for sufficiently small values of e and e’ 



It remains to define F(X) , the field of the initial potential in 

the plane of the variable X, or rather p(e^) in the plane of the 
iX 

variable e . This field is determined by its singularities which are 

IXa iXp 

source vortices in: e and e u as well as in their images with 

respect to the unit circle in the plane of e^ X (fig* IT)* 


The strength of a sink is equal to the strength of a source, and 
may be determined exactly by the condition that the imaginary part of F 
varies by 1 from one vane to the next 


2rtF = (l + i tan a 
(l - i tan a 



iX 

iX 


^iXg-r}^ _ ( 3 _ + i tan p)ln(e iX + e 1 ^" 1 ^) + 

e -iX 3 +T^ _ (1 _ t tan p )ln ( e lX + e " iX 5 +T1 3) 
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The constants tan a and tan 0, put in this form for facilitating 
the ultimate notation, must be such that the complex variable z in the 
plane of the physical variable will be uniform. For this it is necessary 


that the integral 



taken on the contour CDD'C' of the hodograph be 


zero (fig. 15) • Since the function F presents singularities only at A 
and G inside the contour, the latter may be replaced by a loop enclosing 
A and G (fig. 18) . One will note that 


£ ® 5 g « 

where £ and dF j&X are uniform functions of X inside the contour. 

The portions of the integration on the two branches of the loop con- 
necting A and G thus compensate one another, and it suffices to evaluate 
the residues at A and G. 

In the neighborhood of A 

J £ dF ~ £ a ^dF ~ i(l + i tarua)£ A 

since the first logarithm of the expression of F increases by 2jti for 
a contour enclosing A, and the others are uniform. Likewise, in the 
neighborhood of G 


J' £ 3F ~ y'dF i(l +-i tan p)£ G 

The complex variable z will be uniform if 

i(l + i tan a)f; A - i(l +.i tan 0 K g = 0 

C i(q-P) cos P _ Eg 

cos a £ a 

For a deviation by a right angle, the argument of ^ G ^S A is 


a - 0 = 


n 


and 


tan a = -cot 0 = 


Sa 

£a 
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The pitch of the vanes is the variation of Jt, dF for a closed con- 
tour surrounding A or G. This variation has already been calculated above: 

1(1 + 1 tan a)i A = ie 1 * = i(? A + J G ) 


The numerical calculation of the profile, for a real X, offers no 
difficulties. One must especially be careful to choose correctly the 
arbitrary initial values; the holograph method prohibits fixing a priori 
the geometrical values in the physical plane. The angles e and e 1 , 
for instance, determine indirectly the ratio between the chord of the 
profiles and the pitch of the vanes, for it is clear that with short 
profiles considerable velocity differences between infinity upstream and 
the leading edge, on one hand, between the trailing edge and infinity 
downstream, on the other, must be accepted. In compensation, the given 
conditions in the holograph plane permit imposing on the velocity a quite 
continuous variation favorable to a sound flow of the real fluid. 


There exist methods of electrical analogy more powerful and less 
laborious for designing lattices of blades. 

The example above is intended especially for illustrating a study 
of analytic functions, performed in a manner suited to the needs of 
engineers . 


2. ELLIPTIC INTEGRALS GF TEE FIRST KIND AND ELLIPTIC FUNCTIONS 

2.0 Definition 


The elliptic integral of the first kind with real modulus is defined 
by 

u = 


The radical is positive for x = 0. The coefficient k is real and 

smaller than one. Its square k^ is the modulus. The inverse function 
Is denoted by 


L 


cbc 


0 ^(l - x 2 )(l - k 2 x 2 ) 


x = sn(u,k) 

or, more simply, by sn u when k does not vary in the course of the 
calculations. To this function are associated 
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X]_ = cn 
x 2 = 


. U as 


sn^u 

. U = 

l/l - 

k^sn^ 

for 

u = 

0. 

cn u 

and 

dn u 


by integrals 


pi ax 1 

^ 1 JJ[ l _ x.^ ^1 - k^ + k2x i^) 

«--/"* ^ - - 

J 1 ^(1 - * 5 , 2 ) (-1 + k 2 + * 2 ) 


They could therefore be included in the definition of sn u if the 
restrictions fixed for the choice k were lifted; the modulus is only 
under the restriction that it be real. It is more convenient for what 
follows to use the three functions for an equal value of k, smaller than 
one . 

The two f unctions , sn u and dn u are "particular cases of the 
functions t(u) defined in paragraph 1.4 and determine flows around 
doublets in a rectangle. The function cn u likewise defines a flow 
containing doublet s, but with exchange of output between the doublets. 

The three functions are represented by figure 19* The extension of 
the field of cn u by symmetry around a vertical side would complete 
the rectangle along which cn u is real but showing that this function 
is not, like sn u sued dn u, one of the functions t(u) defined in the 
preceding chapter for an arbitrary polygon. Its field contains two doub- 
lets on the perimeter of the rectangle. 


The analytic continuations of the functions sn u, cn u, dn u, 
defined by symmetries around the sides of the rectangles are uniform and 
periodical. The same is true for any rational, function of the three func- 
tions above which is called an elliptic function with real modulus. Such 
a function is not generally real on the contour of a rectangle and does 
no longer define a flow in a rectangle. 

The derivatives of the functions sn u, cn u, dn u result immedi- 
ately from their definitions 
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^•(sn u) = cn u dn u 


— (cn u) = -sn u dn u 
du 


— (dn u) = -k 2 sn u cn u 
du 


In order to normalize the three functions as in the previous chapter 
in linking them to the field of the vortex at the center of the rectangle, 
it suffices to set 


k = sin 0 

sin X 
sin 0 


x = sn u = 


U = / 

Jo 


X 


dX 


\/sin 2 0-. sin 2 X 


/ 

Jo 


dX 


•f 


, 0-X 0+X 0+jt + X 0+Jt-X 

2wsin — - — sin — - — sin sin 


The expressions of dn u and cn u as functions of X are 


dn u = -cos X 

\/sin 2 e - sin 2 X 

cn u = 

sin 0 

The derivative of X(u) is 

— = )4in 2 © - sin 2 X = sin 0 cn u 
du 

The function X(u) is periodic except for a multiple of 2 jt which 
depends on the chosen sections among the vortices. Since the direction 
of the vortices is reversed by symmetry, the periods are twice the sideB 
of the defining rectangle of x(u) , for instance 4K and 4iK* where 
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K* (©) = k(| - ©) 


The function X(u) admits, moreover, the period 2K + 2iK f since 
two successive symmetries reestablish the direction of the vortices. 


2.1 Logarithms of the Elliptic Functions 

The logarithm of an elliptic function, multiplied by a complex coef- 
ficient, is represented by the doubly periodic field of source vortices 
and of sink vortices, corresponding to the poles and to the zeros of the 
elliptic function. The sum of such logarithms is of the same nature but, 
if the coefficients of the logarithms have no simple common measure, the 
intensities of the sources and the circulations of the vortices are varied. 

The logarithms of sn u, cn u, dn u define the fields of figure 21 
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The last field is that of the function sn u, except for one con- 
stant and for one factor, for a ratio of the periods twice that of the 
original function. The conjugate field of the first, with, for instance, 
i cn u dn u/sn u, corresponds to the function dn u for a ratio of the 
periods half that of the original function, in addition, with a change 
of origin. The later calculations will furnish more convenient methods 
for multiplying or dividing the ratio of the periods. 

If the field of the logarithm of an elliptic function is that of 
source vortices, the inverse does not hold true, as the later study will 
show. Even in the case where the representation of a field of source 
vortices is possible by the logarithm of an elliptic function, the finding 
of this function is difficult. The following example which will soon be 
useful shows the origin of the difficulty. 

Let the function 

cn u - cn u, 

F(u) = In =■ 

cn u - cn ug 

be selected for representing a source of the strength 2 ji at u = uj_ 
feeding a sink, or a source of intensity -2rt, at u = ug* 

The field of F(u) contains actually sources of all values of u 
for which cn u = cn u-^ 


u = u^ 


u = u x + 2K + 2iK* 
u = -u-l + 2K + 2iK* 

and at all points deviating from the preceding ones by multiples of 
and of hiK 1 . It is sometimes possible to eliminate the parasitic singu- 
larities introduced in order to conserve only those whose representation 
is desired. 

A particular case is that of a field of source vortices inside a 
rectangle which has already been treated for an arbitrary polygon. The 
logarithm then allows for a constant imaginary part on the contour of 

the rectangle and the transformation z = e^ leads, by analytic con- 
tinuation to the exterior of the circle of radius 1 in the z -plane, to a 
field of source vortices defined in the entire plane. 
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F(u) = 2 A n ln(z(u) - z(u n )] 

where ~ 

2 = dn u + i sin 0 sn u 

If, in particular, the values of the coefficients are real num- 

bers, the function F is the logarithm of an elliptic function. 


2.2 Change of Origin and Theorem of Addition 

The problem of change of origin for an arbitrary elliptic function 
can be reduced to the problem of change of origin of the function X(u), 
that is to say, to the definition of x(u + v) as a function of X(u), 
since the elliptic functions are expressed as functions of X. 

The study may be carried out by conformal transformations, first for 
a real v, then for a purely imaginary v, and finally for a combination. 

It is more direct and more convenient to define the function x(u +_v) 
by its singularities in the u-plane which are vortices with the circula- 
tion *2rt. The conjugate field iX(u + v) is that of sources of intensi- 
ties ±2«. - 


This Is a case where it Is possible to combine the singularities of 
elementary elliptic functions, and one finds 


2ix(u + v) 


= In 


cn u - cn(v - IK') cn(u + iK 1 ) 
cn u - cn(v + iK' ) cn(u - IK' ) 


- cn(v + 2iK' ) 

- cn(v - 2iK' ) 


+ C te 


The singularities of the second term are, in fact, -partly those of 
sources of intensity 2rt at the following points, defined except for 
multiples of 4K and of 4iK' 

v - iK' -v + iK' v + IK' + 2K -v.- iK' + 2K 

v + iK' -v + iK' v - iK' + 2K -v - IK' + 2K 

and partly those of sinks of intensity 2 jt at 


v + iK' 

-v - iK 1 

v - iK' 

+ 2K 

-v + iK' 

+ 2K 

v - iK' 

-v - iK' 

v + iK' 

+ 2K 

-v + iK* 

+ 2K 


still except for multiples of 4K and of 4iK' . 
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Scane of these singularities compensate one another, -whereas the 
effects of the others are additive. The field of the second term is 
definitely that of sources of intensity 4it at 

-v + iK' -v - iK* + 2K 

and of sinks of intensity 4it at 

-v - iK' -v + iK' + 2K 

which are, of course, the singularities of 2ix(u + v). 

If v and u tend successively toward zero, the constant appears 
as zero. 

It remains to link cn(u - iK' ) = -cn(u + iK* ) to X(u) , that is , 
to treat the problem of change of origin initially posed in the parti- 
cular case where v = ±iK* . 

The function i cn(u - iK* ) is real (fig. 24) on the defining 
rectangle of X(u) . It suffices therefore to define it in the plane of 

z = e ix ( u ) where it is represented by doublets at +1 and -1 


i cn(u - iK’ ) = — — + Ci 

z + 1 z - 1 x 

= -cot x(u) + C 2 

The constant C 2 may be specified by 

u = K + iK 1 where cn u = 0 and x(u) = 


It is zero 


i cn(u - iK' ) = -cot X(u) = - 


dn u 


sin 0 sn u 

The formula of the change of origin assumes, consequently, the form 


. v . lain 0snvcnu-idnv^ sin © sn 

2iX(u + v) = In x 

|s 


sin 0snvcnu+idnv sin 0 sn 


u cn v - i dn u 
u cn v + i dn u 


X(u + v) = arc tan Jen v tan X(u)J + arc tan Jen u tan *(v)j 
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From this relationship result the classical addition formulas 

sn u cn v dn v + sn v cn u dn.u 


sn(u + v) = 


2 2 2 
1 - sin 0 sn u sn v 


cn(u + v ) cn u cn v - sn u dn u sn v dn v 
1 - sin^© sn^u sn^v 


dn(u + v) 


dn u*dn v - sin% 6n u_cn u sn v cn v 
1 - sin^0 sn^u s no- 


where the denominator is the square of the common modulus of the four 
terms which appear in the expression of 2ix(u + v) . 

If one has set v = u-j_ and v + u = ~u g, the functions X(u), X(up) 
X(Ug) play symmetrical roles for 


u + u^ + U2 = 0 

It is easy to establish several symmetrical relations between the 
three functions X, but these relations have the disadvantage of not 
defining each of the functions from the two others in a uniform manner. 

The addition formulas are valid for, any arbitrary u and v, but 
for a real u and a purely imaginary v, they reduce the calculation 
of the X function regarding a complex variable to the calculation of 
this function for the real values and the purely imaginary values. It 
is sufficient to change in the formulas v to -iv. 


More conveniently, one may note that the mode of definition of X 
by a symmetrical vortex establishes the relation 

x[k + iv, sin X@j = +> x{v - K’ , cos ©J 

or simplifying the notation 

X(K + iv) = | + X' (v - K' ) 


. , sin [x(K + iv)l cos fx* (v - K' j] dn(v - K' ) 

sn(K + iv) = = — *s 


sin 0 


sin 0 


sin 0 



Y 


MCA m 1^55 


25 


dn(K + iv) = cos £x(K + iv)J = -sinlx* (v - K' )]] = -cos 0 sn(v - K' ) 


cn(K + iv) 


L sn 2 (K + iv) 
I sin^0 


= -i cot 0 cn(v - K' ) 


It is then indicated to replace in the addition formulas u by 
u - K, and v by -K - iv. The functions of u + iv will thus be 
expressed with the aid of functions of u - K and of v - K’ ; the addi- 
tion formulas, in their normal forms, permit transforming these functions. 

The calculation leads to the classical formulas 


X (u + iv) 



i 


arg tanh jcn u tan X* 



snfu + iv) = sn u An 1 v i cn u dn u sn 1 v cn v 

^ * P o o p — — - 

cn* v + sin 0 sn u sn’ v 


cn(u + iv) = 


cn u cn* v - i sn u dn u sn* v dn'v 
cn'^v + sin%> sn^u sn'^v 


dn(u + iv) = 


dn u cn* 


v dn 1 v - 1 sin^Q sn u cn u sn* v 
cn’^v + sin% sn^u sn’^v 


Thus, it is sufficient to calculate x(u) 
functions can be deduced, for the values of © 
and the real values of u between 0 and K. 
formula establishes, moreover, for a change of 
the real period 


from which the other 
between 0 and n/2. 
However, the addition 
origin by a fourth of 


dn(u - K) 


= dn(K - u) 


COS 0 

dn u 


It 


cos {x(K - uj] 
suffices therefore to make 

dn g = \jcos~Q Bn | = ■ 


cos |x(u)] = COB Q 
u vary from 0 


£ 


COS 0- 
2 


eng 


to K/2 where 

\/cos 0 
i 2 cos ^ 
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We remark furthermore that 


cn u 


tanpc(K - u)_ 
tan 0 


and this formula permits calculation of cn u from X(K- u ) , known from 
0 to K, without use of a radical. --- 


The addition formula of X may be written 

x(u + v) - arc tanfes. x(K - v)l + aro tan [tan x(v)tan X(K - u)1 

l_ tan © J |_ tan 0 J 

This expression explains why use of the symmetrical relationships 
between the three functions x(u), x(-v), and x(-u - v) is not very 

convenient. The function X(u + v) is actually expressed with the aid 
of four nonindependent functions. 

One could set up a parallel argument regarding a change of origin 
by a fourth of the imaginary period 


sin X(u + iK* ) sin x(u) = sin 0 

Another classical means for calculating cn u without using a 
radical consists in setting 


cn u = cos cp 


sinh 
tan 0 


sn u s sin cp 


sin X 
sin 0 


dn u = cos 8 cosh i|r = cos X 

2icp = In c , n . q.+Asq.-P. 

cn u - i sn u 


2i{f = In sin 9 dn u + cn u 
sin 0 dn u - cn u 

The fields of the functions <p(u) and >|r(u) are represented by 
figure 25. The functions icp(u) and i{r(u) have fields which present 
the same singularities as ln[sn u] and ln(dn u] except for the ratio 
of the periods and the origin. This could, furnish a second means for 
doubling or dividing the ratio of the periods by two. 
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The derivatives of the functions cp and t|t are 

“£.dnu 

du 


dilr 

— = -sin 0 sn u 
du 


The functions ty, x, 9 therefore integrate, except for one fac- 
tor, the functions sn u, cn u, dn u. 

Changes of origin permit integration of other simple elliptic func- 
tions with the functions X, cp, i|r calculated for u + Kj u + iK 1 J 
u + iK + iK 1 . 


2.3 Gauss and Landen Transformations 


The most direct and geometrically most significant method for 
studying the doubling of the ratio of the periods consists in the study 
of the x(u) function. 

Let us consider two functions x(u,k) and *1 (up,k^ j same 

origin (fig. 26 ). 

The scale of up is chosen in such a manner that the real periods 
of the two functions are the same in the u-plane 

%-! 

The ratios of the periods are 

3-K'l _ i K 1 
Kp 2 K 

In the plane of the variable z = e^, the function Xp is repre- 
sented by two vortices at and -z 0 inside the unit 

circle the circulations of -which are ^2it . It is defined by these singu- 
larities and their images with respect to the circle. 


*1 


i In 


z + Zp 
z - z 0 



+ C te 


The minus, sign was introduced into the logarithm so that the con- 
stant should become zero when Xp is calculated for z = 1 where 

Up = u = 0 
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Xi = 2 arc tan 


2z q sin X 

T 77 ^. 


At 


u = K + IK* , the function X is «/2 and the function X^ is 


rt - 0i 


it - ®n =2 arc tan 


~ 2z ( 
1-2 


9 i 

cot -ir = 


sin 


H-f-)] 


We remark that 


sin IX = sln 0 


iK* . . 

sn = i sin 

2 


sn- XL 


. K* 
cn- -j- 


— = i Jsin 0 


and it appears that 0 t_ is connected with 0 by 

tan^ = sin 0 

whereas the relation between X^ and X is 

©1 X 3 

tan -7T- tan -£■ = sin X 


Deriving this last relationship and calculating the derivatives for 
X^ = X = 0, one establishes that 

K, %' 

K K x 


cos 


2 ®1 


and the transformation is finally 


0 1 

tan -*= tan 


pu 1 , a l n _ e l) j _ sin ^ ui cos2 ^ ^ 


"2" fcemj 1 = wb -75“, tail 

It is possible to deduce- from this the relations of Gauss 


r o ®l p ©i 
l |u^ cos^ tan 


2 x 2 ®1 2 .\ 2 ®1 * 2 ®1 

cos^ ~ + sin* 1 -g- sn* 1 u-^ cos' 1 tan^ 


sn^u-p sin = 
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the use of which is less convenient than that of the trigonometric rela- 
tion between the X f unctions. 

Making use of the Gauss transformation, one can obtain from it many 
others by change of origin and interchange of the axes with the aid of 
the addition formulas. The most frequently used one is the transforma- 
tion of Landen which defines the function Xg (u 2 > kg ) as a function of 

X(u,k) for the same origin but a ratio of the real period doubled 
(fig. 27 ). 

It is clear that the Gauss transformation applied around the point 
Kg + iK' g and for the axes turned by a right angle requires a function 

represented by x(u) to correspond to X 2 except for a displacement 

from the origin by -K/2. This geometric consideration guides the 
calculations . 

First of all, the symmetry of definition of X(u) permits making 
a first change of origin 

X 2 (u 2 , sin © 2 j = | + x£i(u 2 - K 2 - iK' 2 ), sin^| 

The transformation of Gauss defines then the vortex in the rectangle 
of twice the length-width ratio 



sin 

Finally, a second change of origin leads to the function X(u, sin©). 


x(-lU 2 + iK-, - K' 2 )=os 2(^ - %), tan S (| - %) 
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[(-u 2 + iK 2 - K‘ 2 )cos 2 ^ - 2p), tan 2 ^ - ^pjj = + 'g, sin ej - | - 

[(.-*- iK'g)cos 2 (j - + K + iK», tan 2 (i - 


rt 

2 


Setting the two last expressions equal: 


cos 


0 = tan2 (4 “ ~i) = rr 


- sin ©2 


sin ©^ 


2Kg iK f 2 

IT ~ lie 5 


U 2 “ 


cos^ 


u 


M) 


cos- 


>/* 0 2\ 

t “ t; 


It remains only to make the calculations with the aid of the addi- 
tion formulas of X(u) 


tan (f " ^) tan 


X(u 2 , sin © 2 ) 

2 4 


= sin 


[x(u + |, sin e) - |J 

= -dn|u + sin ©J 


Noting that 

dn ^ = \Jcoe 0 = tan^ - 


sn ^ = cos 


fc-S) 


tr 

cn ~ *» sin 


(t-») 


we find that the addition formula of dn u furnishes 


t»[f - ISfcsi. 


sin © 2 ) 


_ dn u - (l - cos 9)sn u cn u 
1 - (l - cos 0)sn 2 u 

The notation of this complicated formula has been simplified. It 
is possible to deduce from it Lander 1 s formulas 

, 2/ 


sn^ug, 


sin 0< 




1 - tan* 


1 + tan‘ 


iiM 
it - ¥) 
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dn(u 2 ^ sin © 2 ) = 


2 tan 


1 + tan £ 


( t-a 


/ . \ L 2 /it ®2\ sn u c] 

“(“2. sin e 2 ) = [l - tan 2 (J - -g-jj 

1 - [l + tanS(| - 


cn(u 2 ; sin © 2 ) 


dn(u. 2 , sin @ 2 ) ~ 


2u 


1 - [l - tan 2 (| - ^)] sn^u 


dn u 


where the elliptic functions sn u, cn u, dn u, written in abbreviated 
form, represent 


sn u = sn(u, sin ©) = sn 


uj, nos 2 (| - , f - - &) 


and so on* 

The transformation of Landen permits approximating the logarithmic 
derivative of dn u to the function sn u for a doubled ratio of the 
periods . 

Jin dn(u> sin 0)J = - (l + cos ©) + 1 + COS 0) , 1 ~ s ^ n .- ^j 

The logarithmic derivative of sn u corresponds to a third type of 
transformation derived from that of Gauss by a change of origin, with 
imaginary complex variable , without interchange of the axes . 

Likewise, the functions and e^ u ^ have fields which can 

be linked to those of. the elliptic functions by means of the ratio of 
the periods double or half. 


2.4 Multiplication or Division of the Ratio of the 
Periods by an Odd Integer 

Due to its analytic continuation by symmetries, the function X(u) 
is represented by a field of vortices of alternate directions. If a 
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vortex on p is retained, in the sense of one of the periods, with p 
"being an odd integer, the field established corresponds again to vor- 
tices of alternate directions, that is, to a function X(u) for a ratio 
of the periods p times as large. The initial function x(u) there- 
fore is the sum of p functions X(u) for ratios of periods p times 
as large and being deduced one from smother "by a displacement of the 
origin. This property does not exist if p is an even number as in the 
case of the transformations of Gauss smd Lamden. 

Let us study first the functions x(u) with the same origin for a 
modification of the imaginary periods (fig. 28), comparable to those of 
the Gauss transformation. 

The function x(u^,k^ ) for which the ratio between the imaginary 

period and the real period is p times as small as the corresponding 
ratio for the function X(u,k) is defined by its singularities. 



Using the addition formula of X(u) and noting that 

4 f • ■) ■ -4 tF- 

one has, since X is odd 

^ 2 

x( u l, k l) = x (u,k) + 2 y (— 1) q arc tan 

q=l 

The constant is zero because, for a convenient determination of the 
arc tangents, all terms are zero for u-^ = u = 0 at the origin. 

The values of u^ and of u are linked, according to figure 28 

^1 u 
K X K 

The constants k^ = sin and kg = sin ©5 are likewise con- 
nected. It suffices to write for u = K the relation between X^ 
and X 


tan x(u,k) 




+ c 


te 
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_ p - 1 


0JL = © + 2 


(-l)*! arc tan — 

c HV' c 


cos 0 


Finally, a differentiation of X]_ with respect to u furnishes 


P - 1 


sin 9-j_ cn^Uj_, sin 0^ v • - [l + tan^x(u, 

sin 0 cn(u, sin ©) K “T ^2^. 


sin 0 


)]=nfe, 


q=l cn^lS^-^ cos © + tan2x(u, sin ©) 

l_p J 


value of KLjJ K may be specified for = u = 0 


Kp _ iK i' _ sin 0 x 
K ^ iKj_ _ sin Q 1 


Y“ (-D q 
/qK* 

q=l cn|JL-, cc 


Let us study now the multiplication of the real period by an odd 
integer p (fig. 29) for functions X^v^ioa) and X(u,k) of the same 

origin. 

The same argument leads to writing 

q=Ejll 

x(u 2 ,k 2 ) = ”y~ (-l)Sc^u + k^ + c te 

p- 1 


That is to say 

X^u 2 ,k 2 ) f X(u,k) + 2 V (-l)^ arc tan 

q=l 


Jtan x(u,k)cn^^, 


The relation between k 2 = sin Qg and k = sin 0 will also be 
specified for u = K and u 2 = pK>> 


P-1 


-€£ = 0+2 


(-1) 4- arc tan tan 0 cn|^-, sin 0 
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and the relation between K 2 and K will result also from the calcu- 
lation of the derivatives for u » U 2 = 0 

sin 82 cn(u 2 , sin 9 g) ^ ^ ^ ^~^_g (.p* ^ * tan2 *) cn (|p sin 

sin 0 cn(u, sin 0 ) K ± + ^2* cn 2^ gin 

P-1 

5 s — 

1+2 ^ (-1)^ cn^p,' sin 

q=l 

The above transformations, combined with those of Gauss and Landen, 
permit linking the functions x(u) for the ratios of the periods which 
have an arbitrary fraction as the quotient. However, the calculations 
are very complicated, and the transformations are especially useful for 
the approximate calculation of the elliptic functions. 


e) 


% _ 1 iK> 2 = sin 0 
K ' p iK 1 sin © 2 


2.5 The Functions cp(u) and i|r(u) 

The functions <p(u) and \|r(u) correspond, like X(u), to alter- 
nating vortices or sources and give reason for use of the same method of 
addition of the singularities . 

The function cp(u) is defined by 


sin <p(u) = sn u cos qp(u) = cn u 

Its addition formula may be deduced from those of the elliptic 
functions. — 

We simplify the notation by designating by the subscript 1 the 
functions of the variable u, by 2 those of the variable v, and by 3 
those of the variable u + v 


sin 


<P 5 = 


= sn* = 


sn^ eng dng + sng cni dn^ 


1 - sin^0 sn n 


snr 


cni enp - snn dnn snp dnp 
cos cp 3 = cn 3 = 1 * t 

1 - sin © sn^ Gia^ 
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2icpj cos cp^ + i sin cp^ ( cn l + i sn l ^ 2 ) ( cn 2 + i sn2 dn]_ ) 

cos cpj - i sin 95 ( cn l “ sn l < ^ n 2)( cn 2 “ 1 sn 2 ^ n l ) 


fcidn 2 l 

+ arc tan dn. 

[cni 2J 

eng 1 


= arc tan Jdrtg tan cp-jJJ + arc tan Jdn-^ tan cp^j 

For completing the symmetry with the formula relating to x(u), one 
may note that 

dn u = i tanjcp(K + iK* - u)j 

The same calculation performed for the function Tjr(u) defined hy 


cosh ijr = u 
cos © 


sinh \|f = tan 0 cn u 


leads to the addition formula 


cn-i cno 

t-. = arg tanh sin © — — — — 
5 dn^ dng 


i - arg tanh^sin © snj_ sngj 
_ 


, arg tanh — & '*** S& - arg tanh — *1 (K I u)timh »g (K I 
|_ sin © J L sin ® 

The method of superposition for multiplication or division of the 
ratio of the periods hy an odd integer may be applied to the func- 
tions cp(u) and ilf(u) 


*(»!,*!) - 9 ^ (-D 9 »(« * ^ *) 




= cp(u > k) + 2 


„ tan cp(u,k)dnfe|l, k* ) 

(-1) arc tan - 


*') 
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(p(u + kj 

= <p(u,k) + 2 YZ ^ 

1=1 

For the function t(u), one must he cautious because this function 
is not zero for u = 0, and in order to avoid introduction of a constant, 
it is convenient to calculate i|r(K - u) which is zero for u = 0. 


q. 


r 


arc tan tan cp(u,k)dn 


eH! 






(-D 


1=- 


P- 1 


_P '1 


*(k - u^kj = YZ *( K “ u + ^r~' k ) 


q . : p - 1 


1 = 


_p-l 


V(K - u,k) + 2 21 (-1)^ arg tanh 

1=1 


tanh y(K - u,k) 


dn 




p- 1 

l =±i 5— 

t( K - ugykg) = (-l)* 1 t(K - u + -^,k) 

^__p-l 


1 = 


P-1 


= t(K - u,k) + 2 ^ (-1) arg tanh 

< 1=1 


tanh t(K - u,k)cn^^,l^ 


*#“) 


2.6 Expansions in Trigonometrical Series 

Since the elliptic functions and their logarithms, the latter except 
for a multiple of 2i«, are doubly periodic functions of u, it is pos- 
sible to eliminate one of the periods by a logarithmic transformation, 
foregoing the symmetry of the roles of the . two periods . _ 
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We choose to eliminate the period 4K by the transformation 

ijtu 

2 si JL « In t t = e~^ 

4K 

When u increases by 4 k, the variable t reassumes the same value 
as the elliptic function F(u) -which is a uniform function of t. When, 
in contrast, u increases by 4 IK' , the variable t is multiplied by 
o -2 Jtlw- 

q c = e and thus is not a uniform function of F since to a single 

value of F there corresponds an infinite number of values of t in 

geometrical progression with ratio q^. 

To the rectangle of the periods of the function x(u), there corre- 
sponds in the t plane a ring-shaped area bounded by the circle of the 
radius 1 and the circle of the radius (fig. 30 ) in which the field 

of the X function is that of four vortices with the circulations ±2it 
in 

ijt ( iK* ) 1 1 3 4 . 3 

e ^ « q in -q in q and in -q 


The neighboring rectangles are transformed into rings reduced or 
increased in the ratio q2, and finally the X function is defined in 
the entire plane by its singularities which are those of alternating vor- 


tices on the real axis of t with real complex variables iq' 
r assumes all integer values, positive or negative. 


r + 


where 


r=« 


ix -21 (- 1 ) 


r t - 
In 


„ . 1 


r=-«> 


t + q 


r+ i 


+ C 


,te 


It will be convenient for what follows to distinguish the positive 
and the negative values of r by writing 


CO 


IX - H (-l) r In 
0 


r + i _-p _ 1 

t - q r 2 q 2 + t + c te 

t + <?*h a - r -i -t 


The constant is zero as it appears when u = X and t = 1. 
Coming back to the variable u 


«-z: 


(-l) r In 


0 




( 


1 + .-«)(!- 


q r + 2 e 1 !!) 
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Since the parameter q is essential 1y smaller than 1, every loga- 
rithm can he expanded in series, for the real values of u. Hie validity 
of the addition of the series is less evident but it is accepted without 
justification 

( r + £X 2S ^f(2stl)l58 -(2s+l)iffi 

iX(u) -2EZ (-D q 35 - • “ 

s=0 r=0 2s + 1 L J 


x(u) = k 


> —a i- 

SO 1 + .j 2 " 1 2 s + 




The function x(u) is thus represented by a trigonometric series, 
the coefficients of which are functions of q, which is itself propor- 
tional to the logarithm of the ratio of the periods. The corresponding 
value of 0 is X(K) 


© 


oo , _ .8 s + i 
4 y~ (- 1 ) - g g 
i=0 2s + 1 l + q 2s+l 


and this alternating series is rapidly convergent. 


For K 1 = K, in particular, the value of q is e -Jt and that of 
0 is jt/k- whence the relation 


(-D 


s 


A b *%)* 


* _ _ 

4 T35 28 + 1 l ~ e -(2s+l )* 


The expansion in trigonometric series of cn u is obtained by 
differentiation 


cn u = 


dX 


2rt 




sin 0 du K sin 0 s _o ^ + 2s+l 


cos (2s + l) 


In particular, the value of K is defined for u = 0 
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K = 


at 


co 


□ s+ 2 


Sin 0 8=0 1 + q : 
2rt 


2s+l 

r 


00 co 


sin 0 s=0 r _Q 


(-D r 4 


>(s + J)(r + £) 


2 

~ _ at 
sin © 

The above series permit calculation of the elliptic functions much 
more easily than by numerical integration of the integrals of the first 
kind. 



It is equally possible to deduce from the first expansion of x 

r 

n 


X = 2 ]T C-1) 3 

r=0 


sin 


arc tan 


2r+i • nu 

d cos — 
2K 


w q. r + 2 sing 

^ — + arc tan — 


itu 


1 + q_ ^ cos 


1 - a r + *cos|| 


= 2 y* (-l) r arc tan 
r=0 


r +i 

sing 


1 - q 


2r+l 


with, in particular. 


9 = 2 2_ (-1) arc tan 
r=0 


2 q 


r + i 


1 - q 


2r+l 


The functions cp(u) and i|r(u) give a basis for comparable calcu- 
lations the principal results of which are 


<p s + 2 
2K 


oo a 

q S 


S =1 1 + q 


2s s 




dn u = — + — 
otr tr 


2* K 1 + q 2s 


COS 


(=f) 



4o 
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00 

♦ = -2 YL 

s=0 


1 



2s+l 2s + 1 

- q 




sn u 


2* 

K sin 0 


00 



s=0 



2s+l 

i - q 



+ 



The value of dn u for u = 0 furnishes an expression of K as 
a function of q. independent of 6 


00 n s 

K = i+aty ■■ 

2 s=l 1 + q. 2s 


- | **11 C-if 1 

^ s=l r=l 


1 + 2 ^ q' 
n=l 


n 


When the real period is much larger than the imaginary period, q. 
is close to 1, and the above expansions do not converge very rapidly. 

One may then go back to a ratio of the periods which is more advan- 
tageous, with the aid of one of the formulas which permit dividing the 
ratio of the periods by an odd integer. One may also utilize the above 
expansions, within their limits of convergence, for imaginary values of 
u of the form u = K + iv (where v is real) up to v = K'/2, then 
continue them by the formula of displacement of a half period. 


3. THETA FUNCTIONS 
3.0 Advantage of the Theta Functions 


The elliptic functions and their logarithms are represented by 
singularities of alternating sensej the necessity for this is clear for 
those which are real, or have a real constant part on the contour of a 
rectangle due to the method of analytic continuation by symmetry. This 
limits the possibilities of calculation by addition of the singularities 
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in the whole of the plane. The difficulty appeared first for the multi- 
plication of the ratio of the periods of the function x(u) by an even 
integer. It was also clear for the functions cp(u) and t(u) and, if 
we had preferred to consider the logarithms of the elliptic functions, 
represented by sources the strength of which is conserved in the course 
of the conformal transformations, the difficulty would have continued to 
exist for the elliptic functions whose singularities are poles, which are 
just as easy to add as the logarithms, but whose rules of conservation 
affect the residues, that is, the properties of the logarithms which inte- 
grate them. 

In order to separate the various singularities which reproduce them- 
selves periodically at the intervals 4 k and 4iK' , without associating 
them automatically with singularities of inverse sense, it is desirable 
to define a function presenting singularities which are all of the same 
sense. 

After separating the essential singularities, we have only to con- 
sider those of the sources and of the poles. These latter are deduced 
frcm the sources , whatever their order may be, by differentiation. It 
suffices therefore to study the sources. 

The theta functions are the analytic functions, the logarithms of 
which are represented by the fields of double infinities of regularly 
spaced sources at the vertices of rectangles. One theta function is 
linked to the study of the elliptic functions with real modulus whose 
ratio of periods Is equal to the ratio of the sides of the rectangle 
defining this theta function. In generalizing according to the rules 
of reduction of the elliptic functions , it is stiff iclent that the pro- 
portion of these ratios be expressed by a simple fraction, a ratio of 
two integers. 

An elliptic function with complex modulus may likewise be 

associated with a theta function corresponding to the field of a double 
infinity of spaced sources with periods the ratio of which ceases to be 
purely imaginary. One sees that, if the projection of one period on 
another is a simple fraction, a ratio of two Integers, the theta function 
will be identical with the theta function associated with elliptic func- 
tions of real modulus, except as will be discussed later. 

The elliptic functions with imaginary modulus have, so far, only 
infrequently been applied in fluid mechanics, and do not come within 
the scope of the present study. 


3.1 Definition of the Theta Functions 

The analytic functions are related to the flow of an incompressible 
fluid, and it is not possible to conceive of an infinite number of sources 



k-2 


NACA m 1^35 


without providing for the disposal of their outputs. The theta functions 
prohibit placing sinks at finite distances, and one must admit the exist- 
ence of an essential singularity at infinity which characterizes the infi- 
nite output of the sources. 

It would be possible to reduce the double infinity of the sources to 
a single one by logarithmic transformation as was done in the preceding 
chapter for the x(u) function, but it is simpler to argue directly in 
the plane of definition of the theta function for illustrating the mul- 
tiplicity of determinations of this function. 

The field of a simple infinity of sources, regularly spaced on a 
straight line (fig. 31 ) > is known and corresponds to the analytic function 


In 



Thus, it is possible to define a theta function by superposition of 
terms of the form 


In 


sin 




cosh (rat J) 


where the denominator ch(rrt was introduced in such a manner that 

the logarithm tends toward zero when r increases indefinitely while 
z and b/a remain constant. 


r=oo 

In 0 = y In 

]r»= — oo 




) 


cosh jtr ^ 
a 


+ C 


ite 


„te 

0 = C X 


Ft 


-4 


r"=- oo cosh rtr ^ 


However, this shows the possibility of defining another theta func- 
tion by interchanging the roles of a and b, and of defining an infinity 
of theta functions — combining the two first ones after various changes of 
origin and adding functions which do not have any singularity at finite 
distance like the circular sine and cosine functions. Hie singularity 
to be retained at infinity is therefore not actually determined. 



NACA TM 1^35 


^3 


It suf flees for the intended applications to choose a single func- 
tion for a given distribution of zeros, and only the first definition 
given by the above equation will be retained. A simple transformation 
yields 

• - ° te x Mr) rj 

The constant is chosen real. 

The theta function admits the period 2a. It is real when z is 
real and -when the real part of z is an odd multiple of a/2. It is 
purely imaginary when the real part of z is an even multiple of a/2. 
Its field is represented by figure 32. 

Since the symmetry of the roles of a and b is agreed upon and 
the scale factor is rather inconvenient, the notations are simplified 

T =f T=^=-ilnq 1 < 1 


COS 2 (rt fr) 

cosh 2 (*r |) 


9(v,q) = 2C 0 sin itv j(l - q 2 * e 25tiv ) (l - q 227 e _23tiV )J 

In order to simplify the natation in the calculations which lead 
frequently to a displacement of the origin by a quarter period, one 
denotes the value of 0 after such a change of origin by 0^. 

Q, (v,q) = e ( v + |,q) = 2C q cos *v XT (l + ( l 2r e 2 * 1 ^) (l + q 2 * e _2n:lv ) 
x \ d. / r=l 


Because of the symmetry of the roles of a and b, masked by the 
choice made for the theta function, the calculations lead likewise to 
utilizing a displacement of the origin by T /2, but the func- 
tions 0 + ^,q^ 021(1 ®i( v + jpS-) 8X6 nol; real - ‘wb.en v is real. 
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The product Q^v + ^,q^e^ rtv is real when v is real. It will be 
designated by ^0 2 (v,q) 811(1 used in the calculations, in preference to 

using 9 1 (v + pq]. 

If the origin is displaced once more- by r/2. 


/ 1 

- i / \ 

00 (7 \ / _ \[ 

xe 2 (v + §,q) 

- 0 o a 2 (l + e- 2 ” lv ) 

+ q 2 " • a,lT )(l + q 21 * e- 2 ” 11 ')) 


= q. 2 e” iriv 0 (v,q) 

•*“ - J - 

Hence, the two formulas 

e^v + |,q) - Ae KlV 0 2 (v,q) 

0 s( v + P* 1 ) = e" 11 ^ 0 1 '(v,q) 

In order to complete the symmetry of these two formulas, one chooses 

x - i-aA. 

Elimination of 0 2 after change of v and v + ^ in the first 
formula furnishes 

0l(v + T,q) = q” 1 e" 2iKV e-^q) 

Finally, a 0 ^ function can be defined by 

0 ? (v,q) = 8 2 (y + |,q) 

The functions In 0 2 and In 0^ differ only by the origin of the 
complex variable and are represented by the same field (fig. 33) . 

Finally, the theta function and the auxiliary functions 0^, 0g, 

and 0^ may be defined by the infinite products 
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e(v,q) = 20 o sin [(l - q 2 * e 2 * 1 ^! - q 237 e" 2jtlV )l 

r=l 

9l(v,q) = 20 o cos itv XT [I 1 + ‘a 2 ** e 2,tlY )(l + q 2 * e " 2n:iv )J 

► 

e 2 (v,q) -0 o q'^TT Rl + q 2 ^ 1 e 2 ” lT ) (l + q 2 ” 1 q- 2 " 1 ^ 

r=0 

9 5 (r,q) - 0 o q'FTX [(l - q 2 ^ 1 e 2 > tlV )(l - q 2 " 1 q’ 2 " 1 ^ 

r =0 

with the formulas of change of origin hy a quarter period, written with- 
out mention of the parameter q 



e (▼ + |) = e 1 (v) 

e i( T + l) 

= -e(v) 



02 ( v + 2 ) = e 3^ 

e 3( T + i) 

= S 2 (t) 

> 

hlcvT 

+ 

> J>^, 

CD 

= iq. 4 e _irtV e 5 (v) 

9j(v + | 

j = iq. ^ 

■ -ijtv . . 
e e(v) 

e i( v + i) 

= <l^e“ l3tY e 2 (r) 

e s( T + 5 


1 

-irfv . . 

e ^(v) 


and for a displacement hy a half period 

0(v + 1) = -0(v) 9l(v + l) = -9^(v) 

e 2 (v + i) = e 2 (v) e 5 (v + l) = e 5 (v) 
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fl(v + t) = -q" 1 e “ 2iltv 9(v) 


0l(v + T ) = q" 1 e~ 2± * V e,(v) 


0 5 (v + r) = -q" 1 e" 21rtV e 5 (v) 

0 2 (v + t) ■ q” 1 e “ 2ijtv 0 2 (v) 

iit~ 

The last formula may be interpreted. The function 0g( v )e T a dmi ts 

the period, t. It is real when the variable v is real or purely imaginary. 
It corresponds to the interchange of the roles of a and b in the initial 
notation and its logarithm is represented by figure 


This remark may be made more specific: 

Each one of the terms of 02 (v, q) admits itself an Infinity of zeros 
and can be represented by an infinite product 


1+4* 


2r+l e 2rtiv 


= ( 1 + 


2r+l\ 

* ) 


.P, 


2v 

(2s + 1) - i(2r + l)£jr 


Hence, the purely formal expression 






TT 


s=s- oo r= 




2v 

(2s + 1) - i(2r + 1)^- 

JX 


This product is not convergent, but an artifice that consists in 
arguing in terms of a derivative of sufficiently high order that the 
coefficients of the factors tend rapidly toward zero when r or s 
increase. It is more convenient to derive a sum by calculating 


a? 


In 





00 oo 

--a z: z 

8 = - oo x* = — oo 


2 

(2s + 1) + i(2r + 1)K*/K 

2v 

(2s + 1) - i(2r + 1)K*/K 




This sum is convergent and shows that the function _0g 

which corresponds to the interchange of a and b in the definition of 
©2 is related to the original function by 
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The integration furnishes 

The periodicity observed above leads to replacing v by v + 

K 

in this relation. The formulas for change of origin then furnish 


A = 


rtK 

K* 


B = 0 


The constant C may be specified for v = 0 



Changes of origin permit deducing from the above formula the group 
of formulas for change of axes 
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3.2 Expression of the Elliptic Functions and of Their 
Logarithms With the Aid of the Theta Functions 

An elliptic function is a rational fraction of sn u, cn u, dn u. 
Let m be the degree of its numerator and n the degree of its denomina- 
tor with respect to the whole of these functions. Except for the case of 
reduction, there exist, in. a rectangle of 4K width and 4iK* height, 

4 m zeros and 4 n poles which can be determined as roots of polynomi- 
als. The numerator, for instance, may be written in the form 


G m( sn u,cn u) + dn u B^,.i(sn u,cn u) 
its roots are among those of 

G^ n 2 (sn u,cn u) - ^1 - k^sn^u) H^^ 2 ( sn u J cn u ) 
This function may be written in the form 


®2m( sn u ) + cn u ^-l^ 811 u ) 

its roots are among those of 

S2m 2 ( sn u ) “ ( x “ sn2u ) h 2m-l 2 ^ sn u ) 

This polynomial of degree 4m in sn u has four m roots and 
every root defines four zeros in a rectangle of the periods. Among these 
16 m zeros, only four m are zeros of the numerator of the elliptic 
function since the method that was followed furnishes the zeros of 


G m (sn u, ± cn u) ± dn u E^^sn u, ± cn u) =0 

The determination of the roots of a polynomial of high degree can 
be difficult unless it is only a matter of numerical calculations. Once 
that determination is made, it is easy to separate the convenient roots. 

If the degree m of the numerator differs from the degree n of 
the denominator, the elliptic function presents four poles of the order 
m - n or four zeros of the order n - m, with the common poles sn u, 
cn u, dn u. Finally, the function presents as many poles as zeros, if 
each one is counted with its order. 
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The logarithm of the elliptic function is represented by sources , 
the intensity of which is the product of 2n and the order of the zero 
or the order of the pole, counted negatively. It appears, intuitively, 
that this logarithm is defined by these singularities like a sum of 
logarithms of theta functions. To demonstrate this, it is just as easy 
to reason directly regarding the elliptic function F(u) -which can then 
be written 


F(u) 



■where a^, . • . , ap are the zeros, b i> . . ., b p the poles after 

eventual reduction of the poles common to the numerator and denominator. 

If multiple poles or zeros exist, the corresponding theta function appears 
■with a degree equal to the order. Writing that F(u) admits the periods 
4K and k±K l , one will obtain, according to the formulas for displacement 
of the origin of the theta functions 


f (u + 4-K) * f(u) 


i|Ur+is) 

f(u +• 4-iK' ) = f (u)e ^ 


with 


r + is = b, + . . . + b- - a, . . . -a^ 

1 Pi P 

The function f(u) has neither pole nor zero at a finite distance. 
Its logarithm g(u) has no singularity at a finite distance and follows 
the simpler law 


g(u + 4-K) = g(u) + 2±j0\ 


g(u +• 4-iK 1 ) = g(u) + ~ (r + is) + 2±x\x 


where and p are integers. 

Neither has the derivative dg/du a singularity at a finite dis- 
tance, and since the derivation of the above formulas shows that it is 
doubly periodic, its continuation to infinity does not introduce any 
singularity at infinity. It can therefore only be a constant. 
Consequently 
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g(u) = Au + B 

where A and B are constants . 

Using this expression in the law for displacement of the periods 
for g(u), one finds 


r + is = -4pK + 4iAK* A = i«A/2K 

Hie sum of the complex variables of the poles minus the sum of the 
complex variables of the zeros is a multiple of the periods which cer- 
tainly proves that it is not possible to represent by an elliptic func- 
tion a function defined by arbitrary singularities as was indicated in 
section 2.1. 


Even though the zeros and the poles are defined except for integer 
multiples of the periods, it is still possible to choose them in such a 
manner that r + is, is zero. In this case A is zero, g(u) and f (u) 
are constants: 



The constant is to be determined for an arbitrary value of the 
variable u. Hie logarithm of the elliptic function F(u) is simulta- 
neously determined. 

In the simple cases where the zeros and poles are obvious, it is 
frequently more convenient not to try to make r + is zero but to 
replace the constant by 


c te 


X e^ 


u 

2K 


Finally, if four zeros or four poles differ only by half periods, 
one will be able to replace the product of the four theta functions by 



Let us apply to the function sn u 
3 (llc) vh08e P oles are zeros of 


which presents the zeros of 
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sn u 


0 


sn Uq q/u - iK* 


&L x e (- 


o - iK* \ 


•C^) 


e[— -j . u-uq 

2S L x e^^ir 


or else, choosing u Q = K, expressing 0 as a function of 0^, and 
noting that sn u, 0, 0^ are real with u 


sn 


u e (3t) -- 6; ft) _ 8 (sg) x e 2 (o) 
’ S 5(4) 8 (t) 8 5® M°> 


Likewise 


cn 


u . !i!il x V!| tou = !s(i} x !^ 

e 3 © e i (o) 8 5 &) ^ 

From the last relation, written for u = K, there results 


ro 

1 CD 

2 

~ °° n 2r+l 

IT 1 - q. 

CD 


r=0 , 2r+l 



L 1 + J 


k 1 = cos © = 


From the relation pertaining to sn u, written for u = K + iK 1 , 
results 


k = sin © = 


^°>i 2 = 4 


»2(o) 


n + 1 21 -) 

r=l 


n(i + s 2 ”- 1 ) 

r=0 


These formulas may be utilized for calculation of 9, jointly with 
the formulas of 2.6. Other, more convenient formulas will be set up 
further on. As in 2.6, when K’ < K, one may interchange the roles of 
the two periods, setting 

_ JtK 

q' = e K’ In q* In q = 

which permits still using the above . formulas for a value of q at most 
equal to e -11 . 

The functions <p, X, i|r may be defined by 
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tan X 


tanh ljr a* 


_ 0 2 (o) 

!(&)_ 

e 5 (o) 


®i(°) 


05(0) 

CVJ 

CD 

CD 

H <» 

/-S 

0 

9 i(& 


e 2 ( °) e 3 (o) e 2 (A) 


3.3 Expansions of the Theta Functions in 
Trigonometric Series 

The logarithms of the 0 g and 0 ^ functions may be expanded in 

trigonometric series like the <p, X, i|r functions, but it is more 
advantageous to directly expand the theta functions in trigonometric 
series. 

Let the function be, for instance, 8 2 which, due to its being 

even and to its periodicity, may be written, for a convenient value of 
the constant C 0 left undetermined in section 3. 


02(v,q) =1 + 2 ^ a s cos 2sjtv = +XZ a s e 2slltv 

S=1 S= —00 

Inserting this expression into the relation 


e 2 ( v + = < T 1 e“ 2ijtv 0 2 (r,q) 


and setting the terms equal, one will find 

a g = a s _i = a ( 2 B-l)+( 2 s- 3 )+...+l . ^ 


r - /— 

e p( v 3<l) =1 + 2 ) q s cos 2smv 


6=1 


and consequently 
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/ \ 03 2 
® 3 (v,qj *= 02V V + t^Q.) = 1+2 ^ (-l) S q s cos 2sjtv 
V ^ 1 8=1 


Likewise 


> 1 (v,q) = sjt e 1 ”^ 0 2 (v + 

, S?+M I .(2»l)ter 


s= -00 


” Uif 

= 2 2 _ <l' ^ cos(2s + l)jtv 

s=0 


e(v,q) = -e-j^v + |,q) 


(-l) 8 ! 


(•*& 


sin(2s + l)jtv 


The value of C Q resulting from these expressions will be retained 
hereafter. Setting, for instance, the two values of 6^_(o) equal: 



The trigonometric series of the theta functions are very convergent 
■and convenient for the calculation of these functions and of the elliptic 
functions. Unfortunately, the theta functions do not satisfy any simple 
rule of addition permitting calculation of their values for complex values 
with the aid of their values for real values. 


The values of k and k' corresponding to a value of q may be 
calculated by the formulas 
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sj~ <j( s+ if 


1 + 2 X a s2 

l 


03(0) 

2 

1 + 2 z:(-i)v 2 

1 

02(°) 

j 

_ 1 + 2 21 a s2 _ 



which result immediately from the expressions .of the functions sn u , 
cn u, dn u, with the aid of the theta functions which have been 
established in the preceding paragraph. 


These formulas are retained for the calculation of k and k' , 
in preference to the formulas of 2.6 and 5-2. 1 

One will note, moreover, that the formula established in 2.6 for 
the calculation of the quarter period K is written 

¥ - [ M °>] 2 


3.4 Elliptic Functions With Given Zeros and Poles 


If the zeros a. 


£Lp and the poles b.^ 


bp of an elliptic 


function are given, obviously, for convenience, satisfying the relation: 


ai *$*.»• Sp — b*^ 1 * * ■+ ^p 

the function is defined by 



where F(Uq) is an arbitrary constant chosen as. the -value of the func- 
tion in a complex variable u 0 different from those of the zeros and 
the poles. The problem is to show that it is possible to express this 
function with the aid of the three functions sn u, cn u, dh u. 
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Let us begin with the simple function admitting two zeros and two 
poles . The relation 


a^ + &2 ~ b-j_ + b 2 


expresses that the logarithm of the function is represented by two sources 
and two sinks at the vertices of a parallelogram. One may also regard it 
as corresponding to the superposition of the field of a sink at the center 
of the parallelogram which is capable of absorbing the output of the two 
sources and of the field of one source feeding the two sinks. Besides, 
if the problem is solved when the center of the parallelogram is at the 
origin, the formulas of addition will permit a change in origin. Thus, 
it suffices to study 

" u l\ 




e2 (s) 

This function has the same singularities as 
choosing as a reference u Q = 21K' 


sn‘ 


2 * 1 


sn 2 g 


and. 


sn' 


2 ^ 


2 U 1 
sn^ -£■ 


on 2 | 


!!(i) 




efi - jaVl + p,) 

\2 k-Kj \2 b-KJ 


which becomes after reduction 
1 1 


sn 


2 U 1 


__2 u 
sn — 

2 


e 2 2 (°) 


2 — o 2 

on the other hand, the formulas of addition establish 


2 1 + dn u. 

The problem is thus solved for the elementary elliptic function 
retained. 

Let us go back to the general problem and suppose that the number 
of the zeros and of the poles p is a power of two. If this is not the 
case, it will still be possible to multiply top and bottom parts by the 
same theta functions. 
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Let us then group the terms two by two: 

- Sr^-’ 

where f (u) is expressed as a function of cn u, sn u, dn u. 
Likewise 


4 - = 4^)4 - = 4 ^) ■ 4 - ** 

Since the sum of the complex variables of the roots is equal to the 
sum of the complex variables of the poles, the method eliminates the 
theta functions, and a rational function of sn u, cn u, dn u remains. 


3.5 Use of Theta Functions 


Theta functions appear useful for the calculation of the elliptic 
functions and for the investigation of elliptic functions with given 
singularities . 


They permit establishing more siaqply the properties of the elliptic 
functions demonstrated in the preceding chapter. All demonstrations 
were, in fact, based on the possibility of addition of singularities, 
and it is much more convenient to isolate these singularities rather 
than search for combinations permitting compensation of the parasitic 
sing ular ities. The formulas for multiplication or division of the ratio 
of the periods by an even or odd integer are written very simply 



and 


gCvyq) 


e(y,i r )e(v H- T,q r ) . . .9 |v + (r - l)Ty^ e ici (r-1) (v 
9l(o,q r ) 0 i( T ,q r ) • • • - l)r,r 2 ] 


Finally, the theta functions permit calculation of the elliptic 
integrals as will be shown in the following chapter. 
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We note that it would, have been possible to study the theta func- 
tions before the elliptic functions, but the method would have been 
particularly for these uniform functions with double periodicity, and 
would have concealed the generalization to the Schwarz functions. 


It-. ZETA FUNCTIONS AND CALCULATION OF THE ELLIPTIC INTEGRALS 
4.0 Definition of the Zeta Functions 


The theta functions were introduced, in the preceding chapter, by 
their logarithms, represented by fields of sources analogous to the 
fields of sources and sinks of the logarithms of the elliptic functions. 
It appeared there convenient to use the theta functions for the calcula- 
tions, in preference to their logarithms; however, it is better to study 
the latter whose derivatives are simpler, in the present chapter. 

A zeta function is the derivative of the logarithm of a theta 
function 


Z(v,q) = ~ ln[e(v,q)] 

and likewise Z Z 2 , Z^ ^ (capital £) for 0^, 02> 

Derivation of the relationships 

In [e ( v + 1, q)] = In £-9 ( v, qf) 

ln[e(v+ r,q[] = In [ 0 ( v, q)J - 2i ic(v + ij 

establishes for the zeta function 

Z( v + 1, q) = Z(v,q) 


Z(v+ T,q) = Z( v, q) - 2i it 

The zeta faction admits the real period 1 and, except for a multiple 
of 2i«, the imaginary period r. 


1 The zeta function thus defined is linked to the function Zn(u,q) of 
Jacobi by Zj(v,q) = 2KZn(2Kv,q). It seemed here more convenient to choose 

the logarithmic derivative of the theta function with respect to its normal 
variable v rather than with respect to the normal variable of the elliptic 
functions u which is sometimes 4Kv. 
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Likewise* the differentiation of the relation 

In [ 0 (▼ + |,q)J = ln|i 05 (v,q)J - xi(v + 

given 

z(v + |,q) = Z 3 (v,q) - ni 

shows that, when v is real ^with Z^(v,q) being real^ , the func- 
tion Z^v + £,qj has a constant imaginary part. 

On the other hand, the function Z(v,q) admits a pole at each zero 
of 0(v, q) . It is purely imaginary when the real part of v is -l/2, 0, 
l/2 since In 0 has, on these straight lines, a constant imaginary part 
when dv is imaginary. The representative field is that of figure 35 * 


4.1 Properties of the Zeta Functions 

In the neighborhood of one of the zeros of 0, for instance v^, 
the function In 0 behaves like ln(v - v^) + C^ e , and Z behaves like 
l/(v - V]_) . It admits therefore a double infinity of simple poles of the 

same sense, in contrast to the elliptic functions which have singularities 
of alternating sense. Just as It has been possible to define the logarithm 
of an elliptic function by superposition of logarithms of theta functions, 
it will be possible to define the elliptic function directly by superposi- 
tion of zeta functions. 

Before examining this essential utilization of the zeta functions, 
it is of interest to study seme of their elementary properties. 

If the zeta function represented in figure 35 is derived graphically, 
the field obtained will present a double pole at the origin, the deriva- 
tive will be real for the real values of v - and for those values of v 
which have t/ 2 as their imaginary part since zeta and dv are real. 

It will also be real on the straight lines corresponding to the real 
parts of v: -l/2, 0, l/2, since both dZ and dv are imaginary on 

these straight lines. The field is that of figure 3 It is defined 
at the interior of rectangles on the contour of which it represents a 
real function and corresponds to an elliptic function. 

This function is, except for one constant and one factor, l/sn^2Kv 
or, in other words, the function p(2Kv) of Weierstrass. We shall never- 
theless prefer to evaluate the derivative of the function Z^(v,q) which. 
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as we have seen, differs from. Z(v, q) only by the origin and one con- 
stant (fig. 37) • 

The field is that of a double pole at the center of a rectangle 
representing, always, except for one coefficient and one factor, the 

function dn^2Kv. 


A| 3 (v,ci)] = A dn 2 (2Kv) + B 
which by integration becomes 

Zj(v,q) = C + Bv+ A f [dn 2 (2Kv)Jdv 

J o 

In the neighborhood of v = t/ 2, the Z^ function behaves like 

— and its derivative like - - — - — On the other hand, dn(2KV) 
v - ^ ( v " i) 


behaves like 


£K(v - §) 


and consequently A = 


The Zj function admits the real period 1, and, writing for simpli- 
fication of the notation 

a/2 


E = r [dn 2 (2Kvj] 2K dv 

J o 


we shall find, writing that Z^(v + l,q) = Z^ ( v, q) , that the constant B 

is -4KE. The origin of the symbol E will appear in if. 4 when the 
elliptic integral of the second kind is studied. 

Continuing the integration 


111 - CT + ** fj ** /J 1 ) + E/K l * dT 2 


and writing that 9 ^ admits equally the real period 1, we shall find 
that C = 0. Consequently 
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and. for v = l/2 


Z 5 (v,<i) = 4k2 ^ [to 2 ( 2 ^) - ij X dv^ 

ta i£r ■ -f 1 P w - i] 


X dv c 


“ iffeS •*/.“* *\T [*>*(■"“) -S'"* 

Besides this property of its derivative, the zeta function has the 
advantage of satisfying a law of addition. The logarithmic differentiation 


CD 

+ 

CD 

1 

1 

CD 

IV) 

rv) 

/“N 

O 

1 ■ 1 1 

e 2 (v)e 2 (v 1 ) ©i 2 (o)e 5 2 (o) 

sn 2 | 

Wv-l) sn 2 (2Kv)J 


furnishes in fact 


/ \ / \ , . 4K cn(2Kv)dn(2Kv)sn 2 (2Kv, ) 

Zfv + v., } + Z(v - v n ) - 2Z(v) -irr— - — ^ ^ 

X i/ \ * sn(2Kv)[sn 2 (2Kv) - sn 2 (2Kv 1 jj 

Interchanging the roles of v and v^ and noting that Z is odd, 
one has I - 7 7J 

. v , v 4K cn { 2Kv n )dn (2KV -, ) sn 2 (2Kv) 

zfv + t_\ - zfv - v.) - 2Z(v 1 ) V 1 / 

\ y ' ' ' ' sn (2Kv^ pn 2 (2Kv 1 ) - sn 2 (2Kv)J 

and by elimination of Zfv - VjJ 

. . cn(2Kv)dn(2Kv)sn^(2Kv-i ) - cn f2Kvi ) dn (2Kvi )sn^ ( 2Kv) 

Zrv +Vi ] = Z(v) + ZU) +2K p — - 

sn( 2Kv) sn f2K v-^ ) jsn 2 ( 2K v) - sn 2 (2KVj_)J 

This is a law of addition for Z(v) analogous tojthat of the 
elliptic functions. 

Applying it to v^ = i, gives 
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Z 5 (V) = Z(v) - 2K 


cn(2Kv)dn(2Kv) 

sn(2Kv) 


and this formula relates simply Zj to Z whereas Z]_ and Zg are 
related still more simply to Z and Zj by a simple change of origin. 

Making a change of origin of t/ 2 in the formula of addition of Z 

Zj (v + V]_) * Zj(t) + Zj(v 1 ') - 2Kk 2 sn(2Kv)sn^2Kv 1 ^sn|2K^V + VjJj 

This formula of addition is simple. In order to be able to apply 
it to the calculation of Zj when v is real and v-^ purely imaginary, 

one must have available a means for calculation of Zj(iv). 

The relation between the values of the zeta function for the imaginary 
values and for the real values of the variable result from, the corresponding 
relationship between the theta functions established in 3*1 


In 


0 2 (iv,e lltT )J 


= +i 


s£ + In _e £ ( + + 


-,te 


which furnishes by differentiation 


^(iv.eH + i Z 2 ( + ‘ ¥ ) 


and by a displacement of origin replacing v by v - i 

For practical purposes, it is more convenient to use the relation 
in the form 

Z 3 [iv ♦ l,. 1 -] = 2p * 1 Zj ^ -*] 


for calculating 
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■where the Zj f unctions appearing in the second term correspond to real 

values of the variables when a and (3 are real. The calculation of 
the Z* function for the complex values of the variable is therefore 
reduced to the calculation of the Z3 function for the real values of 

the variable. The other zeta functions are simply reduced to Zj. 

From the above formulas there results an important relation between 
E, E’ f K, K' 



and according to the formula for differentiation of Zj 

41K 2 [dn 2 (2iKv + K,k) - = ^ + 4K' 2 jdn 2 (2K' v + K , ,k') - |j- 

the laws of addition show that this formula is independent of v and is 
written 


E . El 
K K' 


1 + 


_a__ 

2KK' 


k.2 Decomposition of the Elliptic Functions 

Just as it was possible to decompose the logarithm of an elliptic 
function into logarithms of theta functions and consequently the function 
into a product of theta functions, the poles of an elliptic function may 
be approximated by those of a zeta function for definition of the elliptic 
function by the sum of its singularities. . 

Let us first assume that the elliptic function has only simple poles. 
It may be represented by 
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where all poles bj_, . . . b p are distinct. 

In the neighborhood of a pole b^, the function behaves like 



The same is true in the neighborhood of a pole ■which differs from 
b^ only by multiples of 4 K and of ^iK* . 

On the other hand, in the neighborhood of the same poles 




. e j . C o) ■■ 

L J u - b* Q ^u -_p\ 


Consequently, the function 

** •(■Hr 1 ) 


G(u) = F(u) - C 


A=1 


e*(o) x . 



has no singularity at a finite distance. Its derivative is an elliptic 
function "without singularity, even at infinity, by reason of its double 
periodicity. It can only be a constant and consequently 

Cr(u) = A + Bu 

but G(u) admits the period kK and hence B is zero 
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The elliptic function is thus found to be decomposed into zeta 
functions . 


Let us note besides that F(u) admits the period 4iK* and that, 

when u increases by 4iK, the function z(- — J decreases by 2irt. 

' 4k / 

Consequently 


e' 




^ e .(o) X ... X e(^l^E) 


This indicates that the sum of the residues of an elliptic function 
in its poles is zero and expresses simultaneously for a p that is at 
least equal to 3, a general property of the theta functions which one 
must not apply without observing that, for the chosen form of the 
elliptic function, the sum of the complex variables of the zeros is 
equal to the sum of the complex variables of the poles. 

If same poles of the elliptic function are multiples, the calcula- 
tion is appreciably more complicated. 



qO-1 


Noting that Z(v) behaves in the neighborhood of v = 0 like 

dv a ~ 1 


— — — — , we shall form, for v = U " — \ the elliptic function 
~~ kK 


va 


. (u) _!M + ( ,!)<* rn 9 (^r : ) x • ' • x 


This new elliptic function admits the polq_ u = b^ at most to the 
order a - 1; starting again, we shall progressively reduce the order 
and we shall find the principal part in the neighborhood of each pole. 

The function F(u) differs from the sum of the principal parts in the 
neighborhood of each of its poles only by a constant. The elliptic func- 
tion is thus decomposed into a stun of zeta functions and into derivatives 
of zeta functions. 
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For practical purposes, to calculate the principal part in the 
neighborhood of a multiple pole, one will proceed as for the decomposi- 
tion of rational fractions, replacing every term by a limited expansion. 


4.3 Calculation of the Elliptic Integrals 
An elliptic integral is defined by ^p|y, \/(y - yp)(y - y2) > 

i/(y - y?)(y - *4 y] dy where P is a rational function of the three 


variables . 

The form of an elliptic integral is invariant in the course of an 
inversion or of an inversion followed by a symmetry which can be defined 

fcy 


ax + b 
cx + d 


The critical points Xp, xg, x^, x^ are the images, by inversion 
symmetry, of the critical points yp, yg, y^, y^. 

Let us first assume that the four critical points of the y -plane 
are not on the same circle, that Cp is the circle which passes through 
the points y^, y^, y^, that C 2 is the circle which passes through 

Y4> yp, etc., (fig. 38). Let us draw the circle C which bisects 
C 3 , C 4 , and the circle C’ which bisects Cp, C 2 and let us make an 

inversion with respect to one of the points common to these circles. 

The circles C, C* are transformed into straight lines. The circles Cp, 

C 2 on one hand, C^, on the other, are transformed into circles of 

the same radius and, if the origin is chosen at the intersection of the 
images of C, C*, if the x axis is chosen in the direction of the two 
most distant images, for instance x^, Xg, if the scale is chosen so 

that the distance will be 2 , the critical points are 


xg = -xp = 1 X3 = -x 4 = i 

where k is a number which is generally complex and has a modulus smaller 
than 1 . 

The elliptic integral takes the form 


J Pp x, '/l. - x 2 ,\/l - k^x 2 dx 
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For the method of construction which was followed J k cannot he 
real, hut it can he purely imaginary and in this case the integral is 
reduced to that of an elliptic function hy the transformation 



If k is not purely imaginary, the elliptic integral is reduced to 
that of an elliptic function with imaginary modulus, which exceeds the 
scope of the present study. 

Let us now assume that the critical points in the y-plane are on the 
same circle C, that y^_ and y 2 are consecutive, and that C-j_ and C 2 

are the circles orthogonal to C at y^, y^ on one hand, at y 2 and yjj. 
on the other (fig. 39) • Finally, let Cj he the circle tangent to C^, 

Cg, and to C between y^ and y^. 

If we make an inversion with respect to the point of contact of C 
and C^, the two circles are transformed into two parallel straight lines. 

The .circles and Cg have therefore the same radius, and choosing 

the axis of the abscissas following the straight line transformed from C, 
the origin in the middle between x^, x 2 , an & the scale in such a manner 

that x^, Xg is 2, one obtains 

X 1 = " x 2 = 1 x 4 - ” x 3 = £ > 1 
The elliptic integral thus assumes the form 

J P 2 (x, \jl - x 2 , \jl - k 2 x 2 )dx 

and becomes the integral of an elliptic function by the transformation 

x = sn(u,k) 

If now the four critical points y-^, y 2 , y^, y^ are on a circle C, 
the points y^ and y^ and the points y^ and y 2 (fig. 40) are on 
opposite sides. 

The above construction will this time furnish the following values 
for the images of the critical points 
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and the elliptic integral takes the form 

J' P 3 [x,/(l + x) (l + a, 2 x),/(l - x),(l + a^ijdx 
It is reduced to the previous form by the transformation 

x _ (a g - l)xi g - (a 2 + l) 

" ~ («2 - l ) Xl 2 + (a 2 + l) 

Finally, an elliptic integral is the integral of an elliptic func- 
tion. The latter always has a real modiolus when the four critical points 
of the initial integral are on the same circle or on a straight line. It 
does generally not have a real modulus when the four critical points of 
the initial integral are not on the same circle. To make it be effectively 
a function with real modulus, it is necessary that two circles bisecting 
the four circles passing through three critical points be orthogonal. 

This is true in particular for the integral of section 1.6 defining the 
flow inside an equilateral triangle which was reduced to 


n- f— 32S 

J rn? 

the critical points are 1, j, j2, arL & infinity. A bisecting circle 

passing through J, admits as diameter the degenerated circle con- 

necting 1 with infinity. 


With the elliptic function being decomposed into zeta functions and, 
if multiple poles exist, into derivatives of zeta functions, the integra- 
tion is immediate. To the simple poles correspond logarithms of theta 
functions represented by source vortices. Doublets correspond to the 
double poles, doublets of a high, order to multiple poles. 


A first example corresponds to the addition formula of 

e 3 (v) 


y"|sn( 2 Kv)j J^n 2 K(v + V 3 JJ 2 K dv = — — i 


( 2 Kv x ) 


In 


+ v x ) 


vZ ( Tl ) 


Other examples will be given below. 
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4.4 Elliptic Integral of the Second Kind and 
Flows Around Rectangles 

Before transformation into integrals of elliptic functions, the 
elliptic integrals may be decomposed into integrals of .three kinds . 

The calculations are generally more tedious than in the case of immediate 
transformation into elliptic functions, but the integrals of three kinds 
have a direct advantage for many applications. The integral of the first 
kind is connected with the elliptic functions and was studied in the sec- 
ond chapter. The integrals of the second and of the third kind will be 
studied . 

The elliptic integral of the second kind is defined by 



According to the first chapter, it corresponds to the flow in the 
neighborhood of an open rectangular polygon (fig. 4l) . The calculation 
of this integral is immediate. It is sufficient to put 

x = sn(u,k) 


for finding 


p u 

el(u,k) = / dn 2 (u,k)du 

J n 


The vertices at acute angles toward the flow correspond to u = ±K. 
For u = K 


x = sn(K,k) = 1 


e 



dn 2 (u,k)du = E(k) 


a value already introduced in 4.1, in the study of the derivative of the 
zeta function. 


The vertices at obtuse angles toward the flow correspond to 


u = ±K + iK' 


For u = K + iK' 
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x = sn(K + iK' ) = J 


el( K + iK»,k) = E(k) 


f 

J K 


K+UEC* 


dn 2 (u,k)du 


or, making the change of variable: 


u = K + iK’ + iu* 


dn(K + iK* 
dn 2 (K + iK* 


+ iu* ,k) = -k’ sn(u' ,k' ) 
+ iu’ ,k) = 1 - dn 2 (u' ,k’ ) 


Consequently 

el(K + iK 1 ,k) = E(k) + i[K*(k* ) -E'(k»)] 
The function el(u) Is linked to the function 


Is 

is 


eUu,k) = i Z}(±) 

One finds again that for u = K, -when 
E, and that for u = K + IK* , when 


+ Eli 
K 

Zj = 0, the value of el (u) 
= -in, the value of el(u) 


eZ(K + iK' ,k) = =^- + E + = E + i(K* 

2K K 


E’ ) 


In the plane of the variable u, the field of el(u) differs from 
the field of Z^(u ) only by superposition of a uniform stream which 

displaces the stagnation points to lead them to the zeros of dn u 
(fig. 42). 


Of more interest for the applications is the integral 



2 2 
k Sr 


dx 


whose representative field (fig. 43) Is, according to the first chapter, 
that of a flow around a rectangle. 
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It is easy to normalize (fig. 44 ) , by setting 


k = sin X, 


x = 


sin X 
sin X, 


whence follows the new form of the integral 

/sin 2 X Q - sin 2 X 


r 


sin X r 


dX 


btrt it is more practical to reestablish the connection with the elliptic 
functions by setting x = sn(u,k ) , and this transformation gives to the 
integral the form 



O 

cn du = 



2 2 
dn u - k f 


du 


4r el(u) - 
k d 



“^ z 5 (ic) + d- k ' 2 )? 


2Kk‘ 

One of the vertices corresponds to u = K where Z^ 0 and the 

E - k ,2 K 


value of the integral is 


The middle of the side adjacent to 


the origin corresponds to u = K + iK 1 where = -lit and the value 
of the integral is 


- i* + /l _ k . 2)K±ML m 1 - k ,2 K _ ± E 1 - k 2 K' 


2Kk c 


K 


In the plane of the variable u, the field of the new integral 
again differs from that of only by a uniform stream (fig. 45 ) , this 

time leading the stagnation points to the vertices of the rectangle which 
correspond to the zeros of cn u. 

It is equally interesting to construct the field of sn u in the 
plane of the variable represented by the integral 
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This field is represented by figure 46. 

In a general manner, the field of sn u in the plane of the elliptic 
integral of the second kind, increased by a quantity proportional to u, 
is uniform at infinity and surrounds a rectangular polygon the forms of 
■which vary -with the real proportionality factor. 


4.5 Elliptic Integral of the Third Kind 


The elliptic integral of the third kind is defined by 



dx 

^ 1 - )(l - k^x^) 


In contrast to the two preceding integrals, it depends on two parame- 
ters: k and a. 


Taking into consideration the methods of integration applied above, 
it is ■useful to study the integrals of the third kind only in the cases 
of direct significance for the applications. These latter correspond 
almost always to real values of a which alone are considered here. 

The position of a with respect to the critical points: -l/k, 

-1, 1, l/k determines the appearance of the rectangular polygon, 

corresponding to the real values of a and limiting the representative 
field which, in accordance with the first chapter, is that of a doublet 
of the order l/2. 

Figure 47(a) corresponds to a value of a lying between -1 and 1. 

Figure 47(b) corresponds to a value of a between 1 and l/k. 

Figure 47(c) corresponds to a value of a larger than l/k. 

It is of no interest whatever to normalize here in defining the 

field of a vortex. It is preferable to set 

x = sn(u,k) a = sn(u Q ,k) 

This transformation leads to the integral 


r du 

J 0 sn u - sn Ujj 


which can be easily represented in the plane of the variable u 
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For u 0 between -K and K (fig. 48(a)), the function to be inte- 
grated presents a pole at u = Uq in the neighborhood of which it behaves 
like 



i 

1 


r z v “] 

j - sn u 0 cn u 0 dn Up^u - u 0 ) 

sn 

Uq + (U - Up) 


and consequently the integral behaves like ^ ^ — — ln(u - Uq). It 

is represented by a source. It is real, like its derivative, for 
- K < u < K and for -K + iK* < u < K + iK’ . It is purely imaginary 
along the vertical lines of abscissas -K and +K since its derivative 
Is real and du is imaginary on these straight lines. The derivative 
is zero for u = iK 1 which corresponds to a stagnation point of the 
stream. 

For Uq lying between K and K + IK’ (fig. 48(b)), the same 
analysis defines the field of a vortex since cn u Q is purely imaginary. 

For u lying between -K + iK* and K + iK' (fig. 48(c)), the 
field is again that of a source since cn Up dn Up is' real. 

In order to represent the elliptic function with the aid of zeta 
functions, one must determine all poles in a domain of_ variation of u 
from 4K to 4iK' . These poles are the zeros of sn u - sn u Q : 

u Q u 0 + 2iK' -u Q + 2K -u Q + 2K + 2iK' 

We have already calculated the principal part in the neighborhood 

Of U = Up 



sn u - sn u cn u dn u (u - u | 4K cn u„ dn u_ 
o o o\ o } o o 


Proceeding in the same manner in the neighborhood _of the three 
other poles, one has 


4K cn Uq dn up _ /u - u Q \ /u - uA 
sn u - sn Up \ 4K / $\ 4k J 



The constant can be determined for 


u = Mq + 2K 
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4K cn Uq dn uq ^ 2K cn Uq dn Ug 


sn u - sn u. 


sn Uq 


• <-r) * <-£> - ^ - 

»fif) * K3) * ‘sfe) 


and. the integration is immediate 
dn Uq 


T U cn u p 
J o snu- 


sn xlq 


du = In 


0 


a Me M 

9 iygg/ 9 g\ra; | 

- U o\ Q / u + M a /M 0 A4A 

L \ 4-K / 2 \ it-K / V+K/ 5 W/ 


Q i u - Uq\ q fu - Xlo 

ei -iK ; 9 3 \ttk~ 


u 

4K 


2K cn Uq dn Uq 


\2K/ 3\2K i 


of: 


sn u o 

One reduces the elliptic integral of the third kind to the integrals 


sn u 


.^u - sn^u 


sn 


sn u - sn 


V 


Agreeing to introduce parasitic singularities ■which compensate one 
another, one has 


sn 


u + sn Uq 


sn u 


sn u r 


sn 


u - sn Uq sn^u - sn 2 ^ sn^u - sn^UQ sn^u - sn 2 ^ 


We shall proceed inversely, deducing the two new integrals from, 
that which has just been calculated 


L 


u 2sn u 0 cn Uq 


sn u - sn 


dn u Q n u cn u dn u r u cn Uq da Uq 

du = / 2 2_ du - / tt-vrv <*u 

P / , I sn u 4* sn u n 

J o sn u - sn u 0 u o ° 

- in fell + uf- d ° +± zfej 

/Uq + u\ [_ sn u o 2K \2KJ 

9 V — 2T"j 


Likewise 
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u 2cn u Q dn u Q sn u 


sririi 


sri 




du 



cn u 0 dn u Q r> u cn u_ dn u Q 

" " - ~ ■ du + / — 

sn u - sn Uq J 0 sn u + sn Uq 



but it is more convenient , in order to avoid the reduction of the theta 
functions, to calculate this integral directly by making the change of 
variable u = K - U]_ which leads to a rational integral in dn u^. 

One can likewise observe that the field of the function — 5 — — 

sn^-u - sriniQ 

differs 'from that of the function i only by the ratio of the 

sn u - sn Uq - 

periods, and one can thus reduce the two integrals to one another by the 
formulas of doubling of the periods. 


h .6 Calculation of the Zeta Functions and 


of Various Constants 


The most convenient method of calculation of the zeta functions 
utilizes the expansions in trigonometric series of the theta functions. 

i2 


J-[e(v,q)l II ( - 1)r 5 

Z( T ,q)=^ L J -° 


( r+ i)‘ 


(2r + l)cos(2r + l)jtv 


e(v,q) 


= it 


^(-1)^ l( r + i) sin(2r + l)*v 


and the analogous formulas for Zj_, Z 2 , . 

It can be equally convenient to directly make use of the expansions 
in trigonometric series. It is sufficient to start from expansions of 
the theta functions into infinite products 


In e 5 (v,q) = In 0,(0, q) - 2 In fj (l - q^ +1 ) + XT In (l - + 

J o J r=0 


r=0 


k- 


lull - 


•) 
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Z 3 (v,q) = -2rti 21 

r=0 


2r+l 2iitv 

00 

2r+l -2irtv 

q. e 

+ 2fli T~ 
r=0 

q e 

2r+l 2i:tv 
1 - q e 

n 2r+l . -2iitv 

_1 - q e 


In. expanding, at least for the real values of v, such that 
q 2r+ ^ e*^^ 5TV has a modulus smaller than one: 

_ / \ _ . v 2- v~ (2r+l)s 2isitv , _ . v - (2r+l)s -2isitv 

Z 3 (v,q) = -2*i > 2_ q v ' e + 2rti 2_ 2_ 4 e 

r=0 s=l r=0 s=l 


w _ g 

= lfit — - sin(2sjtv) 

s=l 1 - q 2s 

Likewise 

00 

Z 2 (v,q) = Zj^v + - 4« 21 (- 1 ) 6 £ qS ^ 2s sin(2s«v) 

For Z and Z^_ one must isolate the poles 

V 2 - q2s 

Z(v,q) = jt cot(jtv) + 2 — — r sin(2sitv) 

s=l 1 - q2s 


“ q2s 

Zj_(v, q) = -jt tan(jtv) + 4rt y (-1) sin(2sjtv) 

s=l 1 _ q 28 


The expansions of the derivatives of the zeta functions are deduced 
from the expansions above, and for instance 


dZ 5 

dv 



21 — S - ^o ~ cos(2sjtv) 
s=l 1 - q 


This defines E(q) for v = 0 


K[K - E] _ 


2ji £ 


S 1 


s=l 1 - q‘ 


2s 
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One may also utilize, for the calculation of E, the theta functions 

z - -L 
z 3 “ ej dT 


dZ^ 

dv 


1 d2e 3 1 / d9 5V 

63 dV 2 0* 2 \dv / 


and for v = 0 


8rt 2 

dZ 3 = n=l 
dv 

1 + 2 


(-D^nV 2 


00 _ P 

(- 1 ) q. 


n=l 


K(K - E) = n=l 
2 


(-l) n+i n 2 q n 


2rt 


1 + 2 (-l)V 1 

n=l 


The decomposition of the elliptic functions into zeta functions 
introduces 


dv*- ~ * v = 0 


This quantity may be calculated in various ways, for instance 
starting from the expression of sn u 


9 (ah (o) 

e j(a;) 9 i (o) 


Deriving 


1 

M-k) . 

^'Gk) 6 ^) 

e 2 (°) 

2 K 

_ e 5(ic) 

9 3 2 (i) 

CD 

cT 


and for u = 0 



NACA TM 1^35 


77 


i e'(o) q 2^ 

2 K ej® e^) 


taking into account 



we shall find 


0'(o) =1182(0)92(0)0^(0) 


This formula specifies the expansion of the theta function into an 
infinite product 


QCvjGL) = 92( 0 ) e 2^ 0 ) 6 3^ 0 ) 8 ^ n 1 1 

r=l 


(l - q2 r e 2: *- :rtv ) (l - q^ r e -2 ^ ltv ) 


(X - ^) : 


Finally, the principal constants may be calculated with the aid of 
the following formulas 


{W = e 2(°) 


/k" = '/sin Q = — 


n=l 
00 

ay , 


( n+ 5)" 


00 2 

,n 


0 2W 1 + £Z« 

n=l 


a , , 1 + 2 2 _ (- DV - 

^ = sa 

9 2 <°) . . f=- _2 


USE q“‘ 
21=1 


, E _ 2*2 e 3 "W 
1 * K 


9 ,"(o) 


11=1 


(-l) n+ 1 n 2 g .“ 2 


K 2 9 3 (0) 28 4 (o) 8,(0) 2 (1 + 


00 o\4 CO _ £ 

z; i a£ ) 1 + z: <-i)V 

n=l / n=l 
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It will always be possible to assume in these f ormulas K' ^ K 
and q ^ e“ rt by interchanging eventually the roles of K and K* . 

Moreover, the infinite products of the expressions of the theta 
functions may be reduced to series. 

For instance 


In 


IT a - #*) 

r=l 


,*n(i - q^) 


r=l 


03 00 

-ULi i 2 ™ 


r=l s=l 


= -Hi 


,2s 


s=l 8 1 - q 28 


4.7 Second Integrals of Elliptic Functions 

If one must again integrate the integral of an elliptic function, the 
weighted sum of logarithms of theta functions and of zeta functions or of 
derivatives of these iatter, only the logarithms introduce new functions. 

For the applications and especially for the numerical calculations it 
may be sufficient to define these new functions by series. For instance 


TO 

ln[e 3 (v)] = In 83(0) + y In 

“ "" r=0 


"(l - q 2r+1 e 2iltv )(l - q2r+l e -2ijt) 


(1 _ q 2r+l)‘ 


= m 8,(o) *t: t: 

5 r=0 s=l 8 L J 


In 83(0) + 4 2T 


s=l 


>(l - q 28 ) 


0 

sin sjtv 
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In a general manner, the complicated integrals pertaining to elliptic 
functions may he calculated numerically by use of expansions in series. 

The method does not make evident the singularities of the functions except 
those which are situated on the contour of integration} they are most fre- 
quently led back to the real axis and must be removed from the expansions 
in series. 


5. VARIOUS APPLICATIONS 
5-0 Ranges of Application 


The present paper was written so as to permit understanding and 
facilitate utilization of the studies of fluid mechanics for which we 
introduced elliptic functions. We refer to these studies (l to 5) which 
form a group of applications. 

The elliptic functions solve the problems of flows around obstacles 
the contours of which are schematized by segments of straight lines and of 
stream lines. They furnish the field of nonstationary flow around plane 
vanes in vibration. 

The examples chosen below illustrate the possibilities of use. 


5.1 laminar Profiles for Airplane Wings and Blade Grids 

The problem is to define profiles along which the velocity varies 
very slowly in order to avoid separation of the boundary layer and to 
delay the appearance of shock waves at high speeds. The importance of 
this problem is emphasized in several of our publications and in the work 
done by my coworkers under my direction (6 to 15 ). 

The airplane wing is a particular case of a blade grid for an infinite 
blade pitch. It is therefore sufficient to study the blade grid. 

For the first schematization which will be corrected later on, we 
assume that the velocity is rigorously constant on certain regions of the 
profiles. The pressure is then equally constant for the isentropic poten- 
tial flow, and the contour is that of a free jet in equilibrium with an 
inert space at constant pressure. In order to simplify the calculations, 
the connections between the above regions are assumed to be rectilinear. 

A profile may thus be defined by a dihedron AB, AB' at the leading 
edge, two stream lines BC and B’C', two segments CD, with C'D forming a 
dihedron at the trailing edge (fig. 49). The lower surface of the profile, 
convex near the leading edge and concave in the central region if the cam- 
ber is great, normally admits a point of inflection I. 
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We disregard the possibilities of application to fluids satisfying 
a judiciously chosen law of compressibility which result from a simple 
modification of the integration beginning with the hodograph, and treat 
only the flow of the incompressible fluid (3 and 10 to 15)* 

If F is the complex potential of the flow in the physical plane 
where the complex variable is z, we construct the hodograph £(F) = dz/dF 
or rather the logarithmic hodograph A = ln£. 

Along the stream lines, the intensity of the velocity is constant, 
the modulus of £ and the real part of A are constants. 

Along the segments of straight lines, the direction of the velocity 
is constant, the argument of £ and the imaginary part of A are constants. 

The image of the profile in the plane of the hodograph is therefore a 
rectangular polygon (fig. 50), the general appearance of which can be easily 
constructed from an approximate outline of the flow. 

The field F(A) is that of a source vortex with "the image 0 of the 
infinity upstream discharging into a sink vortex with the image O' of the 
infinity downstream. In the case of a wing profile, this field is reduced 
to the field of a doublet with circulation with the unique image of infin- 
ity. It is not difficult to define as in 1.7, if one knows how to make a 
half plane correspond to the interior of the rectangular polygon by a con- 
formal transformation. ** 

Let us choose to define such a conformal transformation f (A) by the 
field of a doublet at N on C'D (fig. 51)* The reduction of the elliptic 
integrals shows that it is always possible to choose the point N, the ori- 
gin and the scale of the field in such a manner that the function assumes 

the values -1 and 1 at B and B*, l/k and -l/k at C* and C. 

Let us thus assume fj_, f 2 , f 3 to "be the values of f at A, the 

image of the stagnation point, D the image of the trailing edge, I the 

image of the point of inflection at the lower surface of the profile; the 
derivative dA/df then is determined by its singularities except for one 
factor: "• 

dA _ C te JT+Jkf f ~ f 3 

“ " - “ ( f - f i)( f - f a) 

According to figure kl, the constants satisfy the inequalities 

-1 < f ! < 1 - 1 


f 2 < -l/k 
■l/k < f 5 < -1 
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and if one sets 

f = sn u f 1 = sn u^_ f 2 *» sn(u2 + iK’ ) 

the function dX/du is elliptic: 

dX _ (1 + k sn u) (sn u - 13) 

■°° (sn" u I u - sn(u2 "+ iX' )J 

The constant C Q is real and positive because u is real along BA, 
and contained between -K and U]_, -whereas dX/du is real and positive. 

It is convenient here to decompose the function which depends only- 
on sn u. 


_ _ 1 _ d^i _ +Ji + ^1 + ^2 

C Q du sn u - sn u^_ sn u - sn ru2 + iK' ) 

by replacing the parameter f 5 by the parameters a]_ and 02 linked by 
the condition that dX/du should be zero for k sn u = -1 

^ ocq_ 0^2 <x-j^ c Cg sn ii£> 

1 + k sn u^_ 1 + k sn ^12 + iK' ) 1 + k sn u^ 1 + sn ug 

The function X(u) is the sum of integrals of the third kind and of 
a linear function of u. 

Before proceeding with the calculation, it is suitable to remark that 
extensive work would not be justified for a crude representation. A much 
more significant profile may be obtained for a continuous contour of a 
holograph and close to the polygonal contour. It is possible to retain 
two stream lines of the field of a source at A discharging into a sink 
at D (fig. 52 ). The two branches at infinity in A, the image of the 
stagnation point, must be separated from in so that the profile does 
not have an angular point at the leading edge. The two straight lines BA, 
B'A of the polygonal contour must therefore be separated i(n + e + e* ) 
where e and e* designate arbitrary small angles fixing the deviation 
between the holograph of the schematic profile and that of the desired 
continuous profile, and presenting, in particular, a finite curvature at 
the leading edge. 

In the neighborhood of the point A, the derivative dX/du behaves 


like 
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dA 


-a.C 
1 o 


“l C o 


cn u^ dn u-^ u - u-j. 


du sn u - sn uq 
and the integral, except for a regular function, like 


A ~ ‘“ q \ . Q In (u - Un ) 


on u^ dn u^ 


When u exceeds the value uj_ by increasing values, the imaginary 
part of ln(u - u^) decreases by i« and that of A increases hy 
i(n + € + €*); consequently 


“i c ° ■ (*• * H^-) 


Cn 


dn u^ 


Likewise, in the neighborhood of the point D, the derivative dX/du 
behaves like 


d\ 

du 


-q 2 Co k sn U 2 g 2 Cp 

sn u - sn(u + iK 1 ) cn U£ dn u 2 u - u 2 - iK* 


and the integral, except for a regular function, like 

+02 C 0 k sn^U£ 


X Ai 


lnju - Ug - iK*] 


cn U£ dn u 2 

When u exceeds the value u^ + iK* by increasing real variations, 
In Ai - U£ - iK 1 'j increases by +iit, and if A denotes the edge angle of 

the leading edge, X varies by -^1 + ~ - + ^ £l ^ 1A 

_ L t e + e'\A cn u 2 dn u 2 

“ 2 °°-H 1 + — ksn% 

According to the relationship between and established 

before 

C D cn u*^ dn u^_ ^ cn U 2 dn ug 

u = 7 r = C-, 

1 + k sn u*^ * k (1 + sn 112 } sn u 2 


1 + 


€ + € 1 


Hence the expression of the derivative 
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1 + 


£ + e* du 


|A = -kC _ ^ U 1 dn U 1 A 


cn(u 2 + iK‘)dn(u 2 + IK') 


8n 


u - sn u^ 1 11 sn u - sn^U2 + IK') 


Replacing C]_ by its expression and integrating beginning with the 
value Tvg of A at B where u = -K, one obtains 


A - Ab 


1 + 


€ + e> 


= -In 







fet 

- “2 - 

w 

•H 

1 

S' 

1 

4 in 

0 V 

4K / : 

5 \ 4k / 

it 

"Z/!! 

l + U 2 + iK’ Y 

^ f\i + U 2 + iK’ ^ 


L H 

4k / 

2 V 4K ' 


„ 2 K 


u + K 

5k“ 


2 K cn dn u-^ 


1 - k sn u 


sn U]_ 


1 + k sn u-j_ 


j-s)-.© 


A u + K 
it 4 K 


2 K cn u 2 dn u 2 

1 - 

sn u 2 

/u ? \ 

_ 7I d 

_ 2 /V) 

sn 

1 + 

sn U2 

“ All ^ 1 

\ 2 K/ 

3 \ 2 K/ 


The second term is expressed with the aid of theta functions that 
can be calculated for values of the variable which are complex when u 
varies from -K to K by real values. It i6 easy to obtain an expres- 
sion which can be calculated more easily by first coming back to a single 
theta function 


> " u 2 - iK ‘ N 

( /u - u 2 + iK*\ 

V 4K i 

1 V 4k / 

Ju + U2 + 2K + 1K'\ 

+ u 2 + 2K - iK* 

4k / 

A 4K / 


and then decomposing the products of theta functions according to the 
formula of section 3.4 


A 

it 



k + 


k + 


1 + dn(u - U2) 

1 - cn(u - ug) 
1 + dn(u + u 2 ) 
1 + cn(u + U£>) 
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The point B f corresponds to u = K 


V - A B = lrt + K cn U 1 ,il-l sn . u 2 1 A z 


1 + 


e + e 1 


sn uj 


1 + k sc 2 


5 *5 


d) 




cn U 2 dn U 2 1 - sn U 2 
sn Ug 1 + sn U 2 


- A z(^\ - A z C£\ 

2jt \2K/ 2k Z ?V2 K/ 


The real part of Agi - Tvg defines — approximately for the modified 

representation — the logarithm of the ratio of the maximum velocities on 
the upper and lower surface of the profile. 

The calculation of A(u) along the straight lines CA, C x A of the 
hodograph is comparable. The variable u varies by real values from 
-K + iK f to K + IK 1 , and this time it is convenient to make the trans- 
formation of the theta functions containing u^- 

The calculations along the straight lines BC andT B*C T introduce 
the theta functions for complex values of the variable into the two 
logarithms. It suffices to use the formulas of change of given axes at 
the end of 3.1 for arriving at a formula analogous to those used above. .. 

To plot the outline of the modified contour of the hodograph, one 
must use necessarily the theta functions for complex values of the variable, 
but since the values of e and e 1 are small, the calculations converge 
rather rapidly. “ " 

One will begin with defining the image of the field of figure $2 in 
the plane of the variable u (fig. 55). - .r. 

This image g(u) is that of a source at u^ discharging into a 
sink at IK 1 + u 2 with a field which follows the contour of the rectan- 
gle -K, K, K + iK ! , -K + iK 1 . 

The function g(u) is determined by the above singularities. 


g(u) = In 


sn u - sn u-j_ 
sn u - sn (up 4* iK f ) 


= In 


sn u - sn u^ 


sn u * 

k sn U 2 
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If the field of the flow around the profile is not desired, it is 
not necessary to calculate 7\(u) in the entire field of the hodograph. 

The contour suffices for the integration of the profile. However, the 
singularities must he correctly placed. For practical purposes, one must 
start from approximations furnished by a freehand sketch or by electrical 
analogy and then make corrections with the aid of the exact formulas. It 
is, besides, possible to differentiate the latter for corrections in order 
to avoid calculation of the theta f unctions. 

We do not t ake up again the calculation of the integration of the 
profile beginning with the hodograph which is set forth in several of the 
papers cited as references. However, let us note that this calculation 

leads to the use of the functions e^ and e ^ and that the logarithms 
of the theta functions which appear in A are transformed into powers of 
theta functions. 


5.2 Vibrations of a Swallow-Tail Wing 

We have studied the flow around a delta wing in vibration (ref. l6) and 
we intend here to extend the results to a plane, slender swallow-tail 
wing (fig. 5*0* 

We shall not resume the discussion of the approximations which per- 
mit definition of the flow potential by the real paid; of 

where x is the abscissa counted from the vertex, y Q the algebraic 

distance to the axis of symmetry on the wing, r the product of the time 
and the velocity at infinity, W n the normal displacement velocity of a 
point of the wing in the course of the vibration, £ = y + iz a complex 
variable formed with the ordinate y and the height z, measured normal 
to the wing, from a point of the flow, f the complex potential of the 
field of the plane flow of a source of the intensity 1 at a point y 0 on 
the upper surface of a transverse section of the wing with the abscissa x, 
discharging into a sink at the same point on the lower surface. The inte- 
gral is extended to the entire section of the wing which depends on the 
abscissa x. ' 

Between the vertex of the wing and the reentrant point of the 
trailing edge, the calculation is identical with the one made for the 
delta wing, and we plan, essentially, to define the function f(£>yo) 
in the transverse sections which intersect the trailing edge. 
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Since the problem is linear, we shall distinguish the symmetrical 
vibrations for which W n is an even function of y Q from the antisym- 
metrical vibrations for which W n is an odd function of y D . This per- 
mits limiting the integration to a half wing, replacing the func- 
tion f (£,y 0 ) by f (£,y 0 ) + f(^-yo) for "the symmetrical vibrations, 
and by f ( £ , y 0 ) - f ( £ , -y Q ) for the antisymmetrical vibrations. In order 
not to complicate the notation, the f function modified in this manner 
will be represented by the same letter. 

Let us treat first the symmetrical problem and neglect, to begin 
with, the influence of the vortex sheet which is shed by the trailing 
edge. The field of the f functipn is that of two sources at symmetrical 
points of the upper surface discharging into two sinks at the same points, 
but on the lower surface (fig. 55)* The intensities are all equal to 
unity. The condition of Joukowsky imposes the position of the stagnation 
point at B, the image of the trailing edge. 

In accordance with a general method for the problems of flow about 
two segments of straight lines, we shall carry out a conformal transfor- 
mation defined by the field of circulation about the segments AB, A’B' 
(fig. 56) 


? = b sn( Ul ,|) 

where a and b are functions of the abscissa x which are linear in 
the particular case where the leading edge and the trailing edge are both 
straight lines. 

In order to avoid using imaginary values of on the wing, we 

shall make the change 

i(ki' - u) = U! - Ki 

The field £(u) is equally represented by the figure 56 , but one 
must interchange the equipotentials and the stream lines. 


£ = a dn(u,k) 


k 2 = l- 


b£ 

a2 


In the plane of the variable u, the representative field of the 
f function can be graphically constructed (fig. 57) with sufficient 
precision to determine the general appearance . and the singularities. 

The rectangle ABOI corresponds to the quarter of the right lower plane 
of figure 56 where the f function admits as a singularity that of a 
sink at M where u = Uq. 
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analytical extension by symmetries defines the complementary 
singularities, nnd the f function is determined, except for one linear 
function of u, corresponding to a uniform flow parallel to the real 
axis, by 


f =- 


1 

rt 



+ C-jU 


df _ 

1 : 

[> ‘ 

du 

" 2Krt 

\ 2K 


+ Cn 


The Joukowsky condition at the trailing edge imposes that df/du be 
zero for u = K 


2Krt 

and consequently 


Z (- 


. ztl * 

[v 

> 2K/ 

\2 2K J 


+ Ci ■ h Zi (it) + Ci 


Jtf = In 


el 

u + Up\ 

2* / 


'U - Uo\ 

. 2K / 


” K Z l(S0 


+ rtC, 


The constant C Q may be defined by 


Co = i ^ + 


IK* „ /Uq 


i(1k) 


K rtK 

so that f will be zero when £ is infinite and u = iK* . 


However, this purely imaginary value may be neglected. 

The problem is therefore solved by the parametric definition of the 
function f (£) with the aid of f(u) and of £(u). 

Let us now study the influence of the vortex sheet originating at 
the trailing edge. 

Because of the linear character of the- problem, the field induced 
by the vortices may be superimposed on the field calculated above tinder 
the condition that this induced field is not altered by the velocity which 
is normal to the wing. This reservation will be respected if the field 
induced by each pair of symmetrical vortices of the sheet encloses the 
mean position of the wing (fig. 58) and satisfies the Joukowsky condition 
at the trailing edge. 
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The field of these symmetrical vortices may be transposed in the 
plane of the variable u (fig. 59) 'where it appears like the field of 
vortices with a uniform flow, parallel to the real axis, placing the 
stagnation point at the image B of the trailing edge. The corresponding 
function f-j_(u) is determined by its singularities and their images of 

the analytic continuation by symmetry. For a vortex with a circulation 
equal to one 


fl 


_i_ 

2k 


In 




+ D 0 


+ D]_ u 


The stagnation point is at B if df^/du is zero for u = K. 

0 = -i- Z^(—'\ + D-, = — Zn + — — + D-, 

2kK 3\2K) 1 2rtK' - L \2K ' * ^ / 2KK* 1 

where the logarithmic derivatives of the relations established at the 
end of 3.1 have been utilized for the calculation of the function 

for a purely imaginary argument. 

To facilitate the notation, the Z]_ function will be represented 

by 

- V(ir) 

the f-^ function is then defined by 



The constant D 0 can be determined by the condition that f^ should 
be zero for infinite £, that is, for u = iK* . 


Nevertheless, it appears as purely imaginary and may be neglected. 

of the vortex center is 
y D = a dn(K + iK' - ia,k) 

= b sn(a,k') = b sn’a 


The ordinate y Q 
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The equilibrium of pressures on both sides of the vortex sheet 
requires that the potential difference be invariable along a vortex at 
a point -which is displaced at the speed of the flow. 

This variation of potential, equal to the circulation on a contour 
starting from .a point of the sheet and coming back to it, avoiding the 
wing and the sheet, will be the weighted sum of the variations of the 
f and fj_ functions. 

When such a contour does not enclose a vortex, the theta functions 
are uniform and the variation of u is 2K. Consequently, the variations 
of the functions are, respectively 


Af 


=- 2 Z (M 

« H2K / 


Af-, 


ilk 

■ / a- \ 

, ira 


L \2K’ ) 

+ IT 


If the image of the contour in the plane of the variable u encloses 

■i* 


the vortex i In hHr)} it is convenient to add 1 to Af -j 


The potential is the real part of 
'K 


<D = J W n (uo,x,r)f (^u^kjak^n u Q cn Uq dt^ 


r K 

/ 7f^(u,a,k)b cn a dn a da 

Jo ” 


if 7 dy Q is the intensity of the vortex of the magnitude 
dy Q = b cn a dn a da. 


is 


The circulation or variation of the potential on the above contour 


r =-£ 

jt 


o 


W n (^O^X,T 


— ^ak^sn Uq cn Uq dUQ - 


b cn a dn a da + 


piK 


7 b cn a dn a da 
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where a Q is the image of the point — with the ordinate y Q — where the 
vortex sheet intersects the contour. 


The condition of constancy of the 
by the flow is written 


circulation for a point entrained 


0 


dT dr ^ dr 


and, consequently, T has the form r(x - T,y 0 ). 

The two first integrals do not depend on y Q and the latter may he 
written 


f r 

J ?n 


Consequently 


r --d^[ r(x - T ' yo) ] 

This relation utilizes completely the law of the variation of T 
with y D and it suffices to retain the expression of r for y Q = b 

in order to eliminate the last integral 



This integro-differential equation in T which can, besides, be 
transformed into an integral equation by an integration by parts of the 
last integral, is not very manageable in the general case, but for the 
applications usually undertaken, the vibrations are harmonic, and the 
solution T is likewise harmonic if the period of establishment of the 
motion is neglected. One agrees to retain only the real parts 


W, 


n 


W„ 


itOT 


(y> x ) e =W Q [a 


dn u D ,x 


0 


iWT 


r = r„(y)e 1 " (T - x ) 


and the integral equation assumes the form 
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■where 

y = b sn 1 a 

and. a, b, k 2 , K, K' are known functions of- x, the two first ones 
linear for rectilinear leading and trailing edges. 

If the last integral is integrated by parts, it takes the form 

i[[i> «■ a) - r o 0,)] [v(£r) + f - 

if* 6 “' - r oto]c(v(*-) + 

The integrated part is zero because Z^' is zero for a = 0 and 
behaves like g) t in the neighborhood of a - K’ whereas, if there 

exists no singularity of the wing contour rendering the derivative of r 
infinite, the first term behaves like 


ar Q aft m ' a) ( _ } 

db da 


dr, 


db 


— b cn 1 


a dm’ a(K* - a) 


cn* a tends toward zero when a tends toward K* and the product of 
the two terms tends toward zero. The argument shows, moreover, that the 
integral which remains is finite. 

If one notes, on the other hand, that 



= 2K* dn ,2 (a - K' - IK) - 2E> 


- 2k 2 K« 


p 

sn* a 


- 2E 1 


cn ,2 a 



92 


NACA 054 1^35 


the integral equation takes the form 

V*) = -| r**/ w o( u o^ x ) z i(^) ak2sn u 0 cn u o du o - 

i J K [r o (b sn a) - r o (b)J £-2k 2 K’ - 2E' 

This equation may be differentiated in order to furnish 

a function of the values of T 0 from 0 to b. It permits 
the construction of r Q (b) by graphic integration for a form of harmonic 
vibration given by W n (uo,x). For practical purposes, it is more con- 
venient to retain the equation in its first form. 

For the transverse cross section of the abscissa x Q passing through 
the point upstream from the trailing edge and constituting the limit of 
the delta upstream giving rise to an easy calculation, the formula fur- 
nishes r(0) for decomposed elliptic functions and zeta functions cor- 
responding to k = 1. Besides, the quantity r(0) is in that case the 
circulation for a contour passing through the axis of the delta wing. 

Starting from this initial value, the value r(b) may be calculated step 
by step in successive cross sections the abscissa x of which varies 
from x Q up to the abscissa of the wing tips. If an effect F(b) is 

known from 0 to b(x) In a cross section of the abscissa x, the curve 
may be extrapolated by a segment of a straight line of undetermined slope 
dr/db. The integral equation, written in the cross section x + Sx where 
5x is a finite variation is then an equation of the first degree in dr/db. 
The procedure, equivalent to the differentiation of the .integral equation, 
has the advantage of avoiding the calculation of the derivatives of the 
elliptic functions and of the zeta functions which depend on the abscissa x 
by the intermediary of their modulus k^. 

When r(b) and its derivative dT/db are calculated, the potential 
is completely determined by the distribution of intensity of the circula- 
tion 7 linked to dr/db. The calculation of the aerodynamic forces may 
then be carried out as for the delta wing. 

If one has to deal with antisymmetrical vibrations, the calculation 
is parallel and it will be sufficient in this case to perform again the 
determination of the functions f and 

Figure 60 represents the field of the function f(£), an & figure 6l 
its image in the plane of the variable u corresponding^to the function f (u) 
which is defined by its singularities, that is, that of a sink at M and 
that of a vortex with the image I of infinity. 


Ida 


!r (b) 

therefore 
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f = 


-1 In 

jt 






4k 


4k 


- i A In 




■nPHHiF) 


+ C Q + CiU 



sn 

{^4 
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f'u + UqX 


sn( 

) 2 'J 


- Acp(u) + C Q + C-jU 


« n m C L r ia ) to C i T Sa ) , x cn H 

-H*) fc 5 


Uq 


du 


u + u r 


- A dn u + C-i 


du 2rt 


u + Uq 


The constants A and C are to he determined in such a manner that 
stagnation points exist at B and 0 in the plane of the variable u. 

Let us study likewise the function f±(0 characterizing the vortex 
sheet (fig. 62 ) , constructing its image in the plane of the variable u 

(fig. 63). 


f, = In 
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+ B<p(u) + D + D-j_u 


£1 = _L 

du 4rt 
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/ u - K - IK' + 1 


- K - IK’ 


+ ia^ 


sn 


^u - K - iK* + ia^ 




- K - iK 1 - ia 
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^dn'(^ 


- K - iK' - ia N 
2 y 


2k 


sn 


i/u - K - iK 1 -in) 


+ B dn u + Dj_ 
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that 


The constants B and Di are to he determined In such a manner 
dfi/du is zero at B and 0. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 29(a). 


Figure 29(b) 








Figure 33. 


Figure 34 
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Figure 58. 


Figure 59. 







